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Abstract 

In this paper we propose a model of random compositions of maps 
of a cylinder, which in the simplified form is as follows: {0,r) G TxM = 
A and 

/ 9 + r + eu±i{6,r). \ 
t r + ev±i{e,r). R 

where u± and v± are smooth and v± are trigonometric polynomials 
in 6 such that f v±{6,r)d9 = 0 for each r. We study the random 
compositions 

{On, Tn) = O • • • O /^o(6»o, Tq) 

with LOk £ 1} with equal probabilities. We show that under natu¬ 

ral non-degeneracy hypothesis for n ~ the distributions of — vq 
weakly converge to a diffusion process with explicitly computable drift 
and variance. 

In the case of random iteration of the standard maps 

9 + r + ev±i{9). \ 
r + ev±i{9) y’ 

where v± are trigonometric polynomials such that f v±(0) dO = 0 we 
prove a vertical central limit theorem. Namely, for n ~ the distri¬ 
butions of Vn — vq weakly converge to a normal distribution AA(0 ,cj^) 
for cr^ = ^ f (v+(9) — v-(9))^ d9. 

Such random models arise as a restrictions to a Normally Hyper¬ 
bolic Invariant Lamination for a Hamiltonian flow of the generalized 
example of Arnold. We hope that this mechanism of stochasticity 
sheds some light on formation of diffusive behaviour at resonances of 
nearly integrable Hamiltonian systems. 
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^University of Maryland at College Park, vadim.kaloshin@gmail.com 
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1 Introduction 


1.1 Motivation: Arnold diffusion and instabilities 

By Arnold-Louiville theorem a completely integrable Hamiltonian system can 
be written in action-angle coordinates, namely, for action p in an open set 
U C M"" and angle 6 on an n-dimensional torus T” there is a function Ho{p) 
such that equations of motion have the form 

9 = oj{p), p = 0, where u{p) := dpHo{p). 

The phase space is foliated by invariant n-dimensional tori {p = po} with 
either periodic or quasi-periodic motions 9(t) = 9 q + toj{po) (mod 1). There 
are many different examples of integrable systems (see e.g. wikipedia). 

It is natural to consider small Hamiltonian perturbations 

H,{e,p) = Ho{p)+eH,{e,p), eeT\ peu 

where £ is small. The new equations of motion become 

9 = u{p) + edpHi, p = -edeHi. 

In the sixties, Arnold [1] (see also lai) conjectured that for a generic 
analytic perturbation there are orbits {9,p){t) for which the variation of the 
actions is of order one, i.e. ||p(t) — p(0)|| that is bounded from below inde¬ 
pendently of e for all e sufficiently small. 

See O Uni EH EH EZ] about recent progress proving this conjecture for 
convex Hamiltonians. 

1.2 KAM stability 

Obstructions to Arnold diffusion, and to any form of instability in general, are 
widely known, following the works of Kolmogorov, Arnold, and Moser called 
nowadays KAM theory. The fundamental result says that for a properly non¬ 
degenerate Hq and for all sufficiently regular perturbations eHi, the system 
dehned by still has many invariant n-dimensional tori. These tori are 
small deformation of unperturbed tori and measure of the union of these 
invariant tori tends to the full measure as e goes to zero. 

One consequence of KAM theory is that for n = 2 there are no insta¬ 
bilities. Indeed, generic energy surfaces Se = {H^ = E} are 3-dimensional 
manifolds, KAM tori are 2-dimensional. Thus, KAM tori separate surfaces 
Se and prevent orbits from diffusing. 
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1.3 A priori unstable systems 


As an interesting model pQ Arnold proposed to study the following example 

P 

= y + Ho{p,q) + := 

P fl) 

= -^-K y + (cosg - 1)+£ili(p,g, 

^ V ^ V 

harmonic oscillator pendulum 

where q,(p,t G T are angles, p, J G M are actions (see Fig. and Hi = 
(cos q — 1) (cos ip + cos t). 



Figure 1: The rotor times the pendulum 


For £ = 0 the system is a direct product of the harmonic oscillator p = 0 
and the pendulum q = sin q. Instabilities occur when the (p, g)-component 
follows the separatrices Ho{p, q) = 0 and passes near the saddle (p, q) = (0, 0). 
Equations of motion for have a (normally hyperbolic) invariant cylinder 
Ag which is close to Aq = {p = g = 0}. Systems having an invariant 
cylinder with a family of separatrix loops are called an apriori unstable. Since 
they were introduced by Arnold [1], they received a lot of attention both in 
mathematics and physics community see e.g. 0 El [im El Ha EH EH ED]. 

Chirikov [10] and his followers made extensive numerical studies for the 
Arnold example. It indicates that the Pdisplacement behaves randomly, 
where randomness is due to choice of initial conditions near HQ{p,q) = 0. 

More exactly, integration of solutions whose “initial conditions” randomly 
chosen e-close to Ho{p,q) = 0 and integrated over time ~ Ine-time. 
This leads to the J- displacement being of order of one and having some 
distribution. This coined the name for this phenomenon: Arnold diffusion. 

Let e = 0.01 and T = —e“^lne. On Fig. 1.3 we present several his¬ 


tograms plotting displacement of the /-component after time T, 2T, 4T, 8T 
with 6 different groups of initial conditions, and histograms of 10® points. In 
each group we start with a large set of initial conditions close to p = q = 
0, / = J*|^ One of the distinct features is that only one distribution (a) is 

^ These histograms are part of the forthcoming paper of the second author with P. 
Roldan with extensive numerical analysis of dynamics of the Arnold’s example. 
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close symmetric, while in all others have a drift. 



A similar stochastic behaviour was observed numerically in many other 
nearly integrable problems m pg. 370, [ini [2H], see also [SB]). To give 
another illustrative example consider motion of asteroids in the asteroid belt. 
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1.4 Random fluctuations of eccentricity in Kirkwood 
gaps in the asteroid belt 

The asteroid belt is located between orbits of Mars and Jupiter and has 
around one million asteroids of diameter of at least one kilometer. When 
astronoters build a histogram based on orbital period of asteroids there are 
well known gaps in distribution called Kirkwood gaps (see Figure below). 

Asteroid Main-Belt Distribution 

Kirkwood Gaps 


Mean Motion Resonance 

(Asteroid; Jupiter) 3:1 5:2 7:3 2:1 



2.0 2.1 2.2 2.3 2.4 2.5 2.6 2,7 2.8 2.9 3.0 3.1 3.2 3.3 3.4 3.5 
Semi-major Axis (AU) 


These gaps occur when ratio of of an asteroid and of Jupiter is a rational 
with small denominator: 1/3, 2/5, 3/7. This correspond to so called mean 
motion resonances for the three body problem. Wisdom |TT] made a numer¬ 
ical analysis of dynamics at mean motion resonance and observed random 
fluctuations of eccentricity of asteroids. As these fluctuations grow and ec¬ 
centricity reaches a certain critical value an orbit of a hypothetic asteroid 
starts to cross the orbit of Mars. This eventually leads either to a collision 
of the asteroid with Mars or a close encounter. The latter changes the orbit 
so drastically that almost certainly it disappears from the asteroid belt. In 
mi we only managed to prove existence of certain orbits whose eccentricity 
change by 0(1) for the restricted planar three body problem. Outside of 
these resonances one could argue that KAM theory provides stability [32] . 

1.5 Random iteration of cylinder maps 

Consider the time one map of H^, denoted 

Fe ■ (p,g,/,</5) ^ {p,q',l',(p'). 

It turns out that for initial conditions £-cIose to Ho{p,q) = 0, except of 
a hypersurface, one can dehne a return map to an 0(£)-neighborhood of 
{p,q) = 0. Often such a map is called a separatrix map and in the 2- 
dimensional case was introduced by physicists Filonenko-Zaslavskii [18]. In 
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multidimensional setting such a map was defined and studied by Treschev 

[331 lag EH HQ]. 

It turns starting near {p, q) = 0 and iterating until the orbit comes 
back (p, q) = 0 leads to a family of maps of a cylinder 

fe,p,q ■ (I^ ‘dO) (/, (p) G A = M X T 

which are close to integrable. Since at (p, q) = 0 the (p, g)-component has 
a saddle, there is a sensitive dependence on initial condition in (p, q) and 
returns do have some randomness in {p,q). The precise nature of this ran¬ 
domness at the moment is not clear. There are several coexisting behaviours, 
including unstable diffusive, stable quasi-periodic, orbits can stick to KAM 
tori, and which one is dominant is to be understood. May be mechanism of 
capture into resonances [I5] is also relevant in this setting. 

In [22] we construct a normally hyperbolic lamination (NHL) for an open 
class of trigonometric perturbations of the form 

Hi = {cosq — l)P(exp(i(p), exp(A)). 

Constructing unstable orbits along NHL is also discussed in [13|. In general, 
NHL give rise to a skew shift. For example, let E = { — 1,1}^ be the space 
of infinite sequences of — I’s and I’s and a : E —)■ E be the standard shift. 

Consider a skew product of cylinder maps 

F:AxE—)-AxE, F{r, 6] u) = {fu}{r, 9),auj), 

where each fi_j{r,9) is a nearly integrable cylinder maps, in the sense that it 
almost preserves the r -component^ 

The goal of the present paper is to study a wide enough class of skew 
products so that they arise in Arnold’s example with a trigonometric pertur¬ 
bation of the above type (see [22]L 

Now we formalize our model and present the main result. 

1.6 Diffusion processes and infinitesimal generators 

In order to formalise the statement about diffusive behaviour we need to recall 
some basic probabilistic notions. Consider a Brownian motion {Bt, t > 0}. 

A Brownian motion is a properly chosen limit of the standard random 
walk. A generalisation of a Brownian motion is a diffusion process or an 
Ito diffusion. To define it let (hi, E, P) be a probability space. Let R : 

^The reason we switch from the (/, (^)-coordinates on the cylinder to (r, 9) is because 
we perform a coordinate change. 
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[0, +C)o) X —)■ M. It is called an Ito diffusion if it satisfies a stochastic 

differential equation of the form 


dRt = b{Rt) dt + a{Rt) dBt, 


( 2 ) 


where B is an Brownian motion, 6 : M —)■ M and a : M —)■ M are Lipschitz 
functions called the drift and the variance respectively. For a point r G M, 
let Pr denote the law of X given initial data Rq = r, and let E,. denote 
expectation with respect to P^. 

The infinitesimal generator of R is the operator A, which is dehned to 
act on suitable functions / : M —)■ M by 


A fir) = lim 

^ ^ * 4,0 


E,[/(i?,)]-/(r) 

t 


The set of all functions / for which this limit exists at a point r is denoted 
Dffr), while Da denotes the set of all /’s for which the limit exists for all 
r G M. One can show that any compactly-supported function / lies in Da 
and that 


A/(r) = Kr)f+ iX(r)|T. 


(3) 


The distribution of a diffusion process is characterise by the drift b{r) and 
the variance a‘^{r). 


2 The model and statement of the main re¬ 
sult 

Let e > 0 be a small parameter and /> 12 be an integer. Denote by Oi{e) 
a function whose norm is bounded by Ce with C independent of e. 
Similar dehnition applies for a power of e. As before E denotes {0,1}^ and 
a; = (..., Uq, ...) G S. 

Consider two nearly integrable maps: 

/^:TxM —^ TxM 

f , _V ( 9 + r + eu^ff9,r) + Oiie^+^,u) \ .. 

■ V ^ / \ r + ev^fi9,r) + e‘^w^fi9,r) + )' 

for cuo G { — 1,1}, where and are bounded C* functions, 1- 

periodic in 6^, Oiie^^°‘,uj) and Oiie‘^~^'^, u) denote remainders depending on u 
and uniformly bounded in ca, and 0 < a < 1/6. Assume 

max \vii9, r)| < 1, 



where maximum is taken over i = —1,1 and all {6, r) G A, otherwise, renor¬ 
malize e. 

We study random iterations of the maps fi and /_i, such that at each step 
the probability of performing either map is 1/2. Importance of understanding 
iterations of several maps for problems of diffusion is well known (see e.g. 

[211E2]). 

Denote the expected potential and the difference of potentials by 
Eu{9,r) := ^{ui{9,r) +u_i{9,r)), Ev{9,r) := ^{vi{9,r) +v_i{9,r)), 

u{9, r) := ^{ui{9, r) - u_i{9, r)), v{9, r) := ^{vi{9, r) - v_i{9, r)). 

Suppose the following assumptions hold: 

IHOl {zero average) Let for each r G M and z = ±1 we have / Vi{9, r) d9 = 0. 

[HI] ( no common zeroes) For each integer n G Z potentials ni(6',?7,) and 
n_i(6',?7,) have no common zeroes and, equivalently, fi and /_i have no 
hxed points; 

[H21 for each r € R we have V v^{0, r)d0 a(r) ^ 0; 

[H31 The functions Vi{9,r) are trigonometric polynomials in 9, i.e. for some 
positive integer d we have 


Vi{9,r) = ^ v^'^\r) exp27iik9. 

/cGZ, |/c|<c/ 


For u G {—1,1}^ we can rewrite the maps in the following form: 

r f 9 \ _ f 9 + r + eEu{9,r) + euJou{9,r) + Oi{e^~^°‘,u;) \ 

\ r J \ r + eEv{9,r) + euJov{9,r) + e‘^w^g{9,r) + Oi{e‘^~^°',u}) )' 

Let n be positive integer and Uk G { — 1,1}, fc = 0,..., n — 1, be indepen¬ 
dent random variables with Pja;*, = ±1} = 1/2 and = {cuo, • • • 

Given an initial condition (6*o,ro) we denote: 


(0n, Tn) := ^o) = O O • • ■ O f^^{9o, Vq). 


(5) 


[H4] ( no common periodic orbits) Suppose for any rational r = p/q E Q 
with p, q relatively prime, 1 < \q\ < 2d and any 6^ G T 



vi{9 + -,r) 
Q 


^ 0 . 


This prohibits /i and /_i to have common periodic orbits of period \q 
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1H61 ( no degenerate periodic points) Suppose for any rational r = p/q E Q 
with p, q relatively prime, 1 < \q\ < d, the function: 

Evp,g{9,r)= 

fcez 

0<\kq\<d 

has distinct non-degenerate zeroes, where En-^ (r) denotes the j-th Fourier 
coefficient of En(6*, r). 


A straightforward calculation shows that: 

n—1 

On = Oo + nro + s E (Eu{6k,rk) + Ev{0k,rk)) 

k=0 

n—1 

+e'^Uk {u{ek,rk) +n(6'fc,rfc)) C>;(n£^+“) 

k=0 

n—1 n—1 

Tn = rQ + e'^EviOk.rk) + e'^ujkv{6k,rk) + Oi{ne^^°-) 
k=0 k=0 


Even though these maps might not be area-preserving, using normal forms 
we will simplify these maps signihcantly on a large domain of the cylinder. 


Theorem 2.1. Assume that in the notations above conditions [H0-H5] hold. 
Let UeS^ —)-s>0ase—)-0 for some s > 0. Then as e —)■ 0 the distribution 
ofvn^ —vq converges weakly to Rg, where R, is a diffusion process of the form 
(|^, with the drift and the variance 

1 nl 

E2{0,R)d0, (t\R)= / v\0,R)d0. 

Jo 



for some function E 2 , 


defined in (11). 


• In the case that u±i = n±i and they are independent of r we have two 
area-preserving standard maps. In this case the assumptions become 

- [HO] / Vi{6)d6 = 0 for i = ±1; 

- [HI] Vi and n_i have no common zeroes; 

- [H21 V is not identically zero. 

- [H31 the functions Vi are trigonometric polynomials; 

- [H41 the same condition as above without dependence on r; 

- [H51 the same condition as above without dependence on r; 
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A good example is ui{9) = vi{9) = cos 9 and U-i{9) = V-i{9) = sin 6'. 
In this case 

[ E2{9,r)d9 = 0, (^^ = [ v\9)d9 

Jo Jo 

and for n < the distribntion r„ —tq converges to the zero mean vari¬ 
ance normal distribntion, denoted A/'(0, More generally, 

we have the following “vertical central limit theorem”: 

Theorem 2.2. Assume that in the notations above conditions [HO- 
H5] hold. Let UeS^ —>'S>0ase—>-0 for some s > 0. Then as 
e —)■ 0 the distribution ofr^^ —uq converges weakly to a normal random 
variable J/{0, 


Nnmerical experiments of Mockel 


show that no common hxed 


points [HI] (resp. [H4]) is not neccessary for Theorem 2.1 to hold. 
One conld probably replaced by a weaker non-degeneracy condition, 
e.g. that the linearisation of maps f±i at the common hxed point 
(resp. periodic points) are different. 


• In [3l] Sanzin stndies random iterations of the standard maps {9, r) 

{9 + r + X(j){9),r + X4>{9)), where A is chosen randomly from {—1, 0,1} 
and proves the vertical central limit theorem; In [2^ l35] Marco-Sanzin 
present examples of nearly integrable systems having a set of initial 
conditions exhibiting the vertical central limit theorem. 


The condition [H3] that the fnnctions vt are trigonometric polynomials 
in 9 seems rednndant too, however, removing it leads to considerable 
technical difhcnlties (see Section 3.2 and Remark 3.1). In short, for 
pertnrbations by a trigonometric polynomial there are hnitely many 
resonant zones. This hniteness considerably simplihes the analysis. 


• One can replace S = {0,1}^ with = {0,1,..., N — 1}^, consider any 
hnite nnmber of maps of the form (|^ and a transitive Markov chain 
with some transition probabilities. If conditions [H2-H4] are satished 
for the proper averages En of v, then Theorem |2. 1| holds. 


3 Strategy of the proof 

3.1 Strip decomposition 

The main idea of the proof is to divide the cylinder A in strips T x /^, where 
C M, j G E are intervals of size for any 0 < (3 < 1/5. Then we 
will stndy how the random variable — vq behaves in each strip. More 
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precisely, decompose the process rn{(^),n G into infinitely many time 
intervals defined by stopping times 

0 < ni < 77-2 < ..., (7) 


where 


• r„. (ca) is £-close to the bonndary between and for some j G Z 

• is e-close to the other bonndary of either or of and 
rij+i > rii is the smallest integer with this property. 

Since e -C being e-close to the bonndary of with a negligible error 
means jnmp from to the neighbonr interval In what follows for 

brevity we drop dependence of r„(a;)’s on u. 


3.2 Subdivision of the cylinder into domains with dif¬ 
ferent quantitative behaviour 

Fix b > 0 such that 0 < (3—2b < 0.04, small 7 > 0, and Ki := Ki{ui, Vi, M 2 , M 2 ), 
7 = 1 , 2 , depending on functions Uj,Vj, j = 1 , 2 , such that Ki < K 2 and all 
are independent of £. Consider the e^-grid in M. Denote by a segment 
whose end points are in the grid. We distinguish among three types of strips 
Ip. We will have strips of three types as well as transition zones from one to 
another. We define: 


• The Real Rational (RR) case: A strip Ip is called real rational if 
there exists a rational p/q E Ip, with gcd(p, q) = I and 1^1 < d. Clearly, 
there are just finitely many strips of this kind. However, this case is 
the most complicated one and requires a detailed study. 

• The Imaginary Rational (IR) case: A strip Ip is called imaginary 
rational if there exists a rational p/q E Ip, with gcd(p, g) = 1 with 
d < |g| < 


The reason we call these strips are imaginary rational, because the 
leading term of the angular dynamics is a rational rotation, however, 
averaged systems appearing in the previous case are vanishing (see the 
next section). 


We show that the imaginary rational strips occupy an action of 

the cylinder (see Lm. A.l in Sect. 0 . We can show that orbits spend 
small fraction of the total time in these strips and global behaviour 
is determined by behaviours in the complement, which we call totally 
irrational. 
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• The Totally Irrational (TI) case: A strip is called totally irra¬ 
tional if r E ly and |r —p/q\ < with gcd{p,q) = 1 , then |g| > 

In this case, we show that there is a good “ergodization” and 

n—1 / \ n—1 

^UkV (oo + k-j ~ ^ uJkV {9o + fcrg) . 
fc =0 k=0 


Loosely speaking, any rQ E lyC] (M\ Q) can be treated as an irrational. 
These strips cover most of the cylinder and give the dominant contri¬ 
bution to the behaviour of — tq. Eventually it will lead to the desired 
weak convergence to a diffusion process (Theorem 2.1). 


• Transition zones, type I: A zone is a transition zone if there is p/q 
such that gcd(p, q) = f and \q\ < d and it is dehned by the correspond¬ 
ing annuli < \r —p/q\ < K 2 £^^^. 

Analysis in these zones needs to be adapted as “influence” of real res¬ 
onances is strong. 


• Transition zones, type II: A zone is a transition zone if there is p/q 
such that gcd(p, q) = f and 1^1 < d and it is dehned by the correspond¬ 
ing annuli < \'f' — p/q\ < 7 . 

Analysis in these zones requires an adjusted coordinates, otherwise, 
we still study the Totally Irrational and the Imaginary Rational strips 
inside of the type II Transition Zones. 


Remark 3.1. Notice that hniteness of Real Rational strips follows from 
assumption [H3]. If the expected potential is not a trigonometric polynomial 
in 9 this is not true. 


3.3 The normal form 

The hrst step is to hnd a normal form, so that the deterministic part of map 
(|^ is as simple as possible. In short, we shall see that the deterministic 
system in both the Totally Irrational case and the Imaginary rational case 
are a small perturbation of the perfect twist map: 

f 0 \ ^ f 9 + r 
\ r J y r 

On the contrary, in the Real Rational case, the deterministic system will be 
close to a pendulum-like system: 

f 9\ f 9 + r 
\ r J ' ^ y r -|- £E{9 
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for an “averaged” potential E(9,r) (see e.g. Thm. 4.2[ We note that 

this system has the following approximate hrst integral: 


H{9,r) = ^ -e E{s,r)ds, 


( 8 ) 


so that indeed it is close to a pendnlnm-like system. This will lead to differ¬ 
ent qualitative behaviours when considering the random system. Inside the 
Real Rational strips as well as the transition zones we use H as one of the 
coordinates. 

The rigorous statement of these results about the normal forms is given 
Sect. 1^ 


in Theorem 4.2 


3.4 Analsys of the Martingale problem in each kind of 
strip 

The next step is to study the behaviour of the random system respectively 
in Totally Irrational, Imaginary Rational and Real Rational strips, as well as 


in the Transition Zones. This is done in Sections 5T -54 More precisely, we 
use a discrete version of the scheme by Freidlin and Wentzell [20], giving a 
sufficient condition to have weak convergence to a diffusion process as e —)• 0 


in terms of the associated Martingale problem (see Lemma B.l). Now using 
the results proved below we derive the main result — Theorem 2.1 This is 


done in two steps. First, we describe local behaviour in each strip and then 
we combine the information. Fix s > 0. 


By the discrete version of Lemma B.l is sufficient to prove that as £ —?• 0 
any time n < se~‘^ and any (0o,''"o) we have 




n—1 


k=0 


- (^b{rk)firk) + 


- f{ro) 0 , 


(9) 


We define Markov times 0 = no < ui < 77-2 < • • ■ < rim-i < nm < n for 
some random m = m{u) such that each n*, is the stopping time as in ([^. 
Almost surely m{u) is hnite. We decompose the above sum 


m 

E 


k=0 


n-k+i 

2 \ ' „—\e^s 


s=nk 


- { b{rs)f{r,) + 


(T^irs) 


firs) 


and show that it converges to /(ro). 
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3.4.1 A Totally Irrational Strip 

Let the drift and the variance be 


b{r) = / E2{0,r)d9 and cr^(r) = / v‘^{9,r)d9, 

Jo Jo 

where the fnnction E 2 is defined in ( 0 - Let ro be e-close to the bonnd- 
ary of two totally irrational strips and let be stopping of hitting e- 


neighbonrhoods of the adjacent bonndaries. In Lemma 5^ we prove that 

(T‘^{rk) 


n ^-1 

2 \ 


e~ Y, e 

k=0 


A/(rfc) - b{rk)f\rk) + 


-firk) 


-/(ro) = 0 ( 6 ^^+'^), 


for some d > 0 . 


3.4.2 An Imaginary Rational Strip 

Let the drift and the variance be 

9-1 


9-1 


biR{9,r) = - 'Y^ E!2{9 + kr,r) and r) = - ^ -|-fcr, r). 


q 


k=0 


q 


k=0 


Let ro be e-close to the boundary of an imaginary rational strip and let n /3 be 
stopping of hitting e-neighbourhoods of the adjacent boundaries. In Lemma 


5.5 we prove that 


np 


-1 




E 

k=0 


A/(rfc) - ( b{9k,rk)f{rk) + 


cr^{9k,rk) 


fin) 


-f{ro) = 0(e2/^+"), 


As one can see, the limiting process does not take place on a line, since the 
drift and diffusion coefficient depend also on the variable 9. 

Notice that the drift b{9, r) and the variance cr( 6 *, r) both are 6 *-dependent 
functions. In Section 3A we show that time spent in these strips is too small 
to affect the drift and the variance of the limiting process. 


15 











3.4.3 A Real Rational Strip 

Let in the rescaled variable r — p/q = R\/e the drift and the variance be 


q-l 


bnR{e,R) = F{e,R), al^{e,R) = {Rp/qfY.^\e + kR,R), 


k=0 


where F is some function to be dehned in (92). Consider the Real Rational 
case assuming that 


|r — p/q\ < 


|R| < 


that is, that r is close to the “pendulum” domain. In this case, we study the 
process {9qn,Hn) with Hn := Rqn), where Hp/'^{9,R) is an approx¬ 

imate first integral of the deterministic system ([^. In the rescaled variables 
it has the form 

d2 

HP/i{9,R) = - 


where 


VP/^i9)= / Evp^g{s,p/q)ds 


for a properly dehned averaged potential (see Thm. 4.2, (14)). In Lemma 


5.6| we prove that, Hn — Hq converges weakly to a diffusion process Rt with 
t = e^n. 

Notice that the limiting process does not take place on a line. In this 
case it takes place on a graph, similarly as in |2I|. More precisely, consider 
the level sets of the function R). The critical points of the potential 

yp/9(5)) give rise to critical points of the associated Hamiltonian system. 
Moreover, if the critical point is a local minimum of V, then it corresponds to 
a focus of the Hamiltonian system, while if it is a local maximum of Rp/'J, then 
it corresponds to a saddle of the Hamiltonian system. Now, if for every value 
H G E. we identify all the points {9,R) in the same connected component 
of the curve {Hp/‘^{9,R) = iL}, we obtain a graph L (see Figure for an 
example). The interior vertices of this graph represent the saddle points 
of the underlying Hamiltonian system jointly with their separatrices, while 
the exterior vertices represent the focuses of the underlying Hamiltonian 
system. Finally, the edges of the graph represent the domains that have the 
separatrices as boundaries. The process Hn takes places on this graph, and 
so it is a diffusion process on a graph. 


3.4.4 A transition zone 


Finally, in Lemma 5.7 we deal with the Transition Zones of Type I and Type 
H, that is the zones in the Real Rational strips such that Ki < |R| < K 2 e~^^^ 
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Figure 2: (a) A potential and the phase portrait of its corresponding Hamil¬ 
tonian system, (b) The associated graph F. 


and < |i?| < In these strips we study the process {9nq, Hnq) = 

{9nq, H{9qn, Rqn))■ In fhis regime we fix small p > 0 and subdivide each zone 
in sub-strips 


= {He^:\H-Ho\<\R,\sR^-P]. 

We prove that, inside each of one these sub-strips, as £ —)■ 0 the process 
Hn — Ho converges weakly to a diffusion process Rt with t = zero drift 
and the variance: 


9-1 


alz{9,R) = + 


k=0 


3.5 Prom the local diffusion in the rational strips to 
the global diffusion on the line 

In this section we resolve the following problem. In order to combine all 
the previous results, which characterise the local behaviour of the process 
inside of infinitesimally small strips, to determine the global behaviour of 
the process in a (9(1)—strip. First, we prove that the Imaginary Rational 
and Real Rational strips cover a negligibly small part of any (9(1)—strip 
(see Section 0 . Then, one can argue that the process is determined by the 
process in the Totally Irrational strips. 

Notice both the drift biR{9, r) and the variance biR{9, r) at any Imaginary 
Rational strip is given by ^-dependent functions. Our main result (Theorem 


2.1), however, is a diffusion process on a line. To prove that this dependence 
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does not enter into the global diffnsion process we show that the process 
spends inhnitesimal amount of time inside of those strips as follows. 

Lemma 3.2. Let {Ok,rk) = f^^{6Q,rQ), k > 1 be a random orbit defined by 
for some random sequence {ujk}k(^z+- Let n < and 

Tnin) = < k < n : rk belongs to either 

an Imaginary Rational or a Real Rational strip}. 

Then for any p > 0 and e > 0 small enough 

P{TR(n) > pn} < p. 


Proof. Dehne 

bmir) := mui^biR{d,r) and crf^(r’) := ^min^ cr|^( 6 ', r) 

Consider the process Rf^ with the drift biR{r) and the variance cr|^(r). By 
definition this process spends more time in that the process with the drift 
biR{9,r) and the variance a]^{6,r). Moreover, it is a diffusion process on a 
line. Then, using a local time argument, it can be seen that the time spent 
on a given domain is proportional to the size of this domain up to a uniform 
constant. Hence, the time the original process spends in all the Imaginary 
Rational strips is inhnitesimally small compared to the time it spends on 
the Totally Irrational ones. However, the time spent in the Imaginary Ra¬ 
tional strip could be infinite and the argument would not be valid. This 
cannot happen, since if r belongs to an Imaginary Rational strip one has 
that a‘^{6,r) ^ 0. Thus, it is enough to prove that for all imaginary rational 
p/q one has cr^{d,p/q) 7 ^ 0. Indeed, if this is true, then for |r — p/q\ < 
and £ is sufficiently small, one has that a‘^{6,r) 7 ^ 0 by Lemma 

Finally, in the Real Rational case one can use a result from [19] that 
diffusion processes on a graph have well-defined local time. Thus, the time 
spent in all the Real Rational strips is inhnitesimally small compared to the 
time spent in the Totally Irrational ones. Now one can have = 0, but it 
happens just when r = p/q, which follows directly from assumption [H2]. In 
this case, one can see that bRR{6, r) 7 ^ 0 , so that the process is non-degenerate 
and thus the fraction of time spend in the Real Rational strips is less than 
any ahead given fraction. □ 

Lemma 3.3. cr‘^{0,p/q) 7 ^ 0 if p/q is any Imaginary Rational. 

Proof. On the one hand, if d < |g| < 2d this is ensured by hypothesis [H4]. 
On the other hand, if |g| > 2d then <j‘^{9,p/q) = 0 implies: 

v{9 + kp/q,p/q) = t), k = t),-■ ■ ,q — 1. (10) 


3.3 
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Now, since v{6,p/q) is a trigonometric polynomial in 6 of degree d, it can 
have at most 2d zeros, or else be identically equal to zero. The latter case 
cannot occur, since by assumption [H2] we know that 


v'^{9,r)^0 for all rG 


so that v{9,p/q) ^ 0. Thus, v{9,p/q) has at most 2d zeros. Consequently 
equation (10) cannot be satished for all fc = 0, • • ■ , g — 1, since |g| > 2d, so 
that <y‘^{9,p/q) ^ 0. The same argument applies to the Transition Zones. □ 


Combining these facts one can apply the arguments from [20], sect. 8 and 
prove that the limiting diffusion process has the drift 6(r) and the variance 
corresponding to Totally Irrational strips. 


cr^l'r 


3.6 Plan of the rest of the paper 

In Section we state and prove the normal form theorem for the expected 
cylinder map E/. Main difference with a typical normal form is that we need 
to have not only the leading term in e, but also e^-terms. The latter terms 
give information about the drift h{r) (see (0). 


In Section 5.1 we analyse the Totally Irrational case and prove approxi¬ 


mation for the expectation from Section 3.4.1 


In Section 5.2 we analyse the Imaginary Rational case and prove an anal¬ 
ogous formula from Section 3.4.2| 


In Section [5^ we analyse the Real Rational case and prove an analogous 
formula from Section 13.4.31 


In Section 5.4 we study the Transition Zones and prove an analogous 
formula from Section 13.4.41 

In Section]^ we estimate measure of the complement to the Totally Irra¬ 
tional strips and the Transition Zones of type II. 

In Section we present several auxiliary lemmas used in the proof. 


4 The Normal Form Theorem 

In this section we shall prove the Normal Form Theorem, which will allow us 
to deal with the simplest possible deterministic system. To this end, we shall 
enunciate a technical lemma which we will need in the proof of the theorem. 
This is a simplihed version (sufficient for our purposes) of Lemma 3.1 in [5]. 

Lemma 4.1. Let g{9,r) G (T x B), where R C M. Then: 

1. Iflo<l and k^O, \\gkir)e^^^’^^\\cio < |fc|'°“^||^||c'o■ 
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2. Let gk{r) be some functions that satisfy < M\k\ " ^ for all 

a < Iq and some M > 0. Then: 


Y1 9k{ry 


27rik0 


kez 

0<k<d 


< cM, 


C'-o 


for some constant c depending on Iq. 

Let TZ be the finite set of resonances of this map, namely, 

n = {p/q e Q : gcd(p, q) = 1, |g| < d}. 

Denote by 0{e) a function whose C°-norm is bounded by Ce for some C 
independent of e. 

Define 

E2{9,r) =Ev{9,r)deSi{9,r)+Ew{9,r), b{r) = [ E2{9,r)d9, (11) 


where Si is a certain generating function defined in (|21|-|22|). 
m ^ 

Ef 


Theorem 4.2. Consider the expected map Ef of the map (|^ 

9 \ _ f 9 + r + eEu{9, r) + 

r J ' ^ \ r + eEv{9, r) + e^Ew{9, r) + 

Assume that the functions Eu{9,r), Ev{9,r) and Ew{9,r) are C\ I > 3 and 
7 > 0 small. Then there exists K > 0 independent of e and a canonical 
change of variables: 

$:TxM ^ TxM, 

{9,h) ^ {9,r), 

such that: 

• If \r — p/q\ > 'J for all p/q eH, then: 

$-1 oE/o<|>(0,h) = 

9 + r + eEu{9,r) + eEi{9,r) + \ (12) 

r + e^E2{9,r) + Oi{e^+<^) )’ 

where Ei and E 2 are some C'‘~^ functions. Moreover, E 2 verifies II-E 2 IIC 0 < 
K and it verifies: 




6 (r) := / E2{9,r)d9 = 

Jo 

drEv{9,r)dgSi{9,r) — dQSi{9,r) (Eu{9,r) —Ev{9,r 


(13) 


d9. 


In particular, b{r) satisfies ||&||co < K and in the area-preserving case 
(when Eu{9, r) = En(6', r) = En(6 *)), h{r) = 0. 
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( 14 ) 


jf |f: — p/q\ < 27 for a given p/q eTZ, then: 

oE/o 4>(0,f) = 

/ 6 + r + £ Eu{ 0 ,p/q) -Ev{6,p/q) + Evp^q{6,p/q) + E3{6) +C)z(e^+“) 
\ f + £Evp^q{ 9 , f) + £‘^E4{0, r) + C)/(e^+“) 

where E.Vp^q is the function defined as: 


E,Vp^q{9,r) = Ew 


kf~\^27vike 


and E 3 is the ^ function: 


Esie) 

where TZp^q = {k E Z : 
Moreover, E 4 is a 


= - Z 

k 0 lp,q 


i{Ev^y{p/q) 

27ik 


27rik0 


k 0 , \k\ < d, kp/q G Z}. 
function verifying ||-E' 4 ||co < K. 


Also, $ is C^-close to the identity. More precisely, there exists a constant M 
independent of e such that: 


$ - Id||c 2 < Me. 


(15) 


Proof. For each p/q G TZ we will perform a different change. Since the 
procedure is the same for all p/q E TZ, from now on we fix p/q E TZ. The 
procedure for the rest is analogous. 

We will consider the canonical change dehned implicitly by a given gen¬ 
erating function S{9,r) = 9r + eSi{9,r), that is: 

r = dgS{9,r) = r + edeSi{9,r) 

9 = dfSi9,r) =9 + edfSil9,r). 


We shall start by writing explicitly the hrst orders of the e-series of o 
E/ o $. If {9,r) = ^{9,r) is the change given by the generating function S, 
then one has: 


$( 0 ,r) = 

9 - edrSii9,r) + e^dedfS4i9,r)dfS4{9,r) + Oie^d^gdfSiid,Sif\\co) \ ( 16 ) 
h + edgS4i9,r) - e^d^sSii9,r)dfS4{9,r) + (!l(e3||a|5i(9,5i)2||co) J ' 


Its inverse is given by: 


$- 1 ( 0 ,r) = 

9 + edrS,{9,r)-eW,S4{9,r)deS4{9,r) + 0{e^dlS4{deS4nco) \ (17) 
r - edgS,{ 9 ,r) + e‘^dgdrS 4 { 9 ,r)dgS 4 { 9 ,r) + OiemdgdjS^ideS^nco) J ' 
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Now, first we compute E/ o $( 6 ',r). One can see that: 


where: 


E/o$(0, h) 


9 + reAi + £^^2 + £^^3 \ 

r + eBi + e‘^B 2 + e^B^ J ’ 


Ai = Ku{9,r) — dfSi{9,r) + deSi{9,r) 

A 2 = —d0Eu(9, f)drSi(9, r) + drEu{6, f)dgSi{6, f) + dodfSiiO, f)dfSi(9, f) 
-d^eSi{e,f)drSi{e,r), 

As = O(||a2a,^i(a,^i)2||co) +0(11^3^1(5,^^ 

+0(||En||ci||5,5i||cii|5,5i||co) + 0(||En||c2(||5e5ii|co + ||5,5ii|co)2), 

and: 


B 2 

Bs 


Ev{d,r) + d 0 Si{d,r), 

—deEv{9, f)dfSi{9, f) + drEv{9, f)deSi{9, r) 



+0(||Eni|c2(||5e^i||co + ||5,^i||co)'). 


Then, using (p(^ one can see that: 


$-3oE/o$(0, r) 


9 + r + eAi + £^^2 
r + eBi + e‘^B 2 + e^Bs 




(18) 


(19) 


where: 


Ai = Ai + dfSi{9 + f,f), 

A 2 = ^2 + £^3 + O(||56i5,S'iy4i||0o) + 0(||5?S'i-Bi||co) 

+Oi\\dlS,dgS,\\co), 

and: 

B, = B,-deS,{9 + r,r) 

B 2 = B 2 - dgSi{9 + r,r)Ai - dfdeSi{9 + r,r)Bi 

+dgdfSi{9 + r,r)deSi{9 + f,f), ( 20 ) 

Bs = 53 + 0(||5,52,Si(5,^i)2||co) 

AO{\\d‘lSi{A2 + £^3)11(10 + \\dedfSiB2\\co) 

+0{^dlSiA\\^co + ll^g^rS'i^iSillco + 115^5?S'l-B^llco) 
+0(||525,^iAi5,5i||co + \\ded^,SiB,dgS,\\co) 
+O(||5,5,^i52^i/li||co + \\{dgdfSiAB,\\co). 
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Now that we know the terms of order e and of o E/ o $, we shall 
proceed to hnd a suitable Si{9,r) such that these terms are as simple as 
possible. More precisely, we want to simplify the second component of (19). 
Ideally we would like that Bi = 0. Namely, we want to solve the following 
equation whenever it is possible: 


deSi {9, r) + Ev{9, f) — deSi {9 + r, r) = 0. 


One can easily hnd a solution of this equation by solving the corresponding 
equation for the Fourier coefficients. To that aim, we write Si and En in 
their Fourier series: 

( 21 ) 

Ev{9,r)= Er;^(r)e“. 

kez 

0<\k\<d 

It is obvious that ior k > d and /c = 0 we can take S'f (r) = 0. For 0 < k < d 
we obtain the following homological equation for S'^(r): 

27iikS'l{k) (1 - + En^(r) = 0. (22) 

Clearly, this equation cannot be solved if = 1, i.e. if kr G Z. We note 
that there exists a constant L, independent of e, L < d~^, such that for all 
0 < A: < d, if r 7 ^ p/g satishes: 

0 < \r — p/q\ < L 


then kr ^ Z. Thus, restricting ourselves to the domain |r — p/q\ < L, we 
have that if kp/q ^ Z equation (22) always has a solution, and if kp/q G Z 
this equation has a solution except at r = p/q. Moreover, in the case that 
the solution exists, it is equal to: 


^f(r') = 


lEv^ir) 


27rk (1 — 


We will modify this solution slightly to make it well dehned also at r = p/q. 
To this end, let us consider a function p{x) such that: 


p{x) = 

and 0 < p{x) < 1 if a; G (1, 2). Then we dehne: 

' _ ^Tvikf 


[ 1 

if 

X < 1, 

1 0 

if 

IV 


k'kif) = h 


2'Kk'y 
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and take: 


5f(f) = 


iEn^(r)(l — /ifc(r)) 


27rfc(l — e 


2TTikf\ 


(23) 


We note that this fnnction is well dehned since the nnmerator is identically 
zero in a neighbonrhood of r = p/g, the nniqne zero of the denominator (if 
it is a zero indeed, that is, if fc e More precisely, we claim that: 


1 if /c e and |r — p/g| < 7 / 2 , 

hfc(h) = ■^ 0 if /c G TZp^q and |r — p/q\ > 87 , 

0 if /c ^ TZp^q. 


(24) 


Indeed if A: G Tlp,q there exists a constant M independent of r and £ snch 
that: 


— |r — p/g|(l — Mir — p/g|) < 
7 


1 — e 


27Tikr 


27ik'y 


< —|r — p/g|(l + Mir — p/g|). 
7 


Then, on the one hand, if k E T^p^q and |r —p/q\ < 7/2 we have: 


^27vikr 


2 TTk'-f 


1 M 

<2 + ^7<1. 


for £ snfhciently small, and thns /rfc(r) = 1. On the other hand, if |r —p/q\ > 
87 then: 


~,27vikr 


27ik'j 


> 8 - 9 M 7 > 2 , 


for £ snfhciently small, and thns /ifc(r) = 0. Finally, if fc ^ '^p,q then: 


-,2'Kikf 


2TTk'y 


> M 7 > 2 


for 7 snfhciently small and then we also have pki'^) = 0 . 


Now we proceed to check that the hrst order terms of (19) take the form 


(12) if |r—p/g| > 87 and (14) if |r—p/g| < 7 . On the one hand, by dehnitions 


(28) of the coefficients S'J^(r) and ( | 20 | ) of Bi, we have: 

Bi= 

0<|A:|<d 


Then, recalling (24) we obtain: 

0 


if \r — p/q\ >7 


Bi = < ^ Er^(h)e™ = Enp,g(0,h) if |h-p/g|< 7 / 2 . (25) 
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where we have used the dehnition (15) of ¥,Vp^q{6,r). On the other hand, 
from the definition ([2^ of S'f (r) one can check that: 


i(En^)'(r)(l - /ifc(h)) + ^ 2 ^ike 


-dfSi{9,r) + dfSi{9 + r, r) 
---d0Si{9+ f,f) - ^ 

0<|fc|<d 


2 'Kk 


Recalling dehnitions (pd| of Ai and (20) of i?i, this implies that: 


Ai = Eu{9, r) — Ev{9, r) + Bi 


- E 

0<|fc|<d 


i(En^)Xh)(l - /ifc(h)) + iEv^{f)ii'Af) ^2nike 
2 Tik 


(26) 

(27) 


Then we use (25) and (24) again, noting that /^^(r) = 0 in both regions 
\r—p/q\ > 87 and |r—p/g| <7/2. Moreover, we note that for Ir—p/gl <7/2. 

Evp^q{9,f) = Evp^q{9,p/q) + 0(7), 

(En'=)'(h) = (En'=)'(p/g) + 0(7). 

i{EvA'{ 


Dehne 


0<|/c|<d 


2 'Kik 6 


2 'Kk 


(28) 


Then the same holds for En(6',r) and En(0,r): recalling definition (15) of 
E3, equation (26) yields: 


^1 = 


En(0, f) — Ev{d, r) + Ei{9, f) 


if |r-p/g|>37, 


i^E{e,p/q)+Evp^q{e) + Ey,{e) + 0{e^lA if \r-p/q\<l/2, 


(29) 


where Eu{9^p/q)— Ev{9^p/q) = AE(6',p/g). In conclusion, by ( [M| ) and ( |^ 
we obtain that the hrst order terms of (17) coincide with the hrst order terms 
of (12) and (14) in each region. 

For the terms we rename B 2 in the following way: 


E2{9,r^ -^ 2 1 {|r— p/ij|>37} ) 

E4{9,r) -®2|{|r— p/ij|<7/2}- 


(30) 

(31) 


Now we shall see that E 2 verifies (13). To avoid long notation, in the following 
we do not write explicitly that expressions Ai, Bi, Ai and Bi are restricted to 
the region {\r—p/q\ > 87}. We note that since in this region we have i?i = 0 
by (25), recalling the dehnition (20) of Bi it is clear that Bi = d0Si{9 + r,r). 
Hence, from dehnition (20) of i?2 it is straightforward to see that: 

B2 = B2-d^0SA9 + k,r)A4. (32) 
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Now we recall that Ai = Ai + dfSi {6 + r,r). Then, using (29), for \r—p/q\ > 
3^^/® we obtain straightforwardly: 


Ai = Eu{9, r) — Eu{9, r) + Ei{9, r) — dfSi{9 + r, r). 


(33) 


Using this and the dehnition (18) of B 2 in expression (32) one obtains: 

^ 2 |{|r-p/g|> 3 ei/ 6 } = -deEv{9, f)dfSi{9, f) + drEv{9, f)d0Si{9, f) (34) 

-dlS^{9,f)dfS^{9,f) 

—d^Si {9 + r, r) [Em( 6 ', f) — Em( 0, r) 

+Ei ( 0 , r) - dfSi (0 + f, r)]. 


Show that this expression has the claimed average (13). On the one hand, 
it is clear that 9gS'i(0, f)dfSi{9, f) and + r, f)dfSii9 + r, f) have the 

same average, so: 

[ -d^eSi{9, f)drSi{9, r) + d^gSi{9 + h, r)drSi{9 + r, r)d9 = 0. 

Jo 

On the other hand, writing explicitly the zeroth Fourier coefficient of the 
product, one can see that: 


-d 0 Ev{ 9 , r)dfSi{9, r) — dgSi{9 + r, r)Ei{9, r)d9 = 0. 


'0 


Thus, recalling (30) and using these two facts in equation (34) we obtain: 

[ E2{9,r)d9 = 




drEv{9,r)doSi{9,r) — dQSi{9,r) (Eu{9,r) — Ev{9,r) 


d9, 


so that (13) is proved. 

We note that, from the dehnition (23) of the Fourier coefficients of Si, 
it is clear that Si is with respect to r. Since it just has a hnite number 
of nonzero coefficients, it is analytic with respect to 9. Then, from the 
dehnitions (30) of E 2 and (31) of E 4 and the expression (20) of B 2 , it is clear 
that both E 2 and E 4 are Finally we shall bound the C°-norms of the 

functions E 2 , b{r) and F ^4 and also the error terms. To that aim, hrst let us 
bound the norms of Si and its derivatives. We will use Lemma 4T and 
proceed similarly as in |S]. We note that: 

1. If Pkir) 7 ^ 1 we have |1 — and thus: 

1 


o27rikr 


< M-V 1^1 
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2. Then, using that ||/ o^||ci < C\\f\^^^J\ci (l + 11^11^;), we get that: 


^ _ ^27rikr 

for some constant M. 


< M 7 


-(Z+l)/5|^|-«-l 


C' 


3. Using the rule for the norm of the composition again and the fact that 
||/i||ci is bounded independently of e, we get: 

for some constant M, and the same bound is obtained for || 1 — pfc(r) ||(^i. 

Using items 2 and 3 above and the fact that ||En^||ci are bounded, we get 
that: 


df 


1 — (f) 

27r/c(l - 


< 


"1 E 


1 


C° ai+a2=a 


|1 - fJ-k{r)\\c° 


1 


1 — e 


2'Kikr 


C“2 




Then, by item 2 of Lemma [4.1[ we obtain: 

11^1 lie* < M7-('+^)/®. 

One can also see that ||9f5'i||c! < M||S'i||c!+i and ||96iS'i||ci < M||S'i||ci. In 
general, one has: 

liax^^illc* < (35) 


Now, recalling definitions (30) of E 2 and (31) of U 4 , and using either 
expression (34) or simply (20), bound (35) yields that there exists some K 
such that: 

||^2||co IIU4IICO < 


To bound h{r) we use again (35) and recall that ||En||co < K, ||5rEn||co < K, 
||Em||co < K. Then from its definition (13) it is clear that: 

ll&llco < K^- 

Similarly, one can easily bound the error terms in the equation for r: 

^34 = 0(£2+-), (36) 

and the error terms for the equation of 9\ 

52^2 = (37) 
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This finishes the proof for the normal forms (12) and (14) (in the latter case, 
we have to take into acconnt the extra error term of order cansed 

by the 7 -error term in (29)). 

To prove (15), we jnst need to recall (16) and nse (35). Then one obtains: 


||$-Idi|c2 <M'5 ||^i||c3. 


□ 

From now on we will consider that our deterministic system is in the 
normal form, and drop tildes. 


5 Analysis of the Martingale problem in the 
strips of each type 

5.1 The Totally Irrational case 

First of all we note that in this case, as in the IR case, after performing the 
change to normal form, the n-th iteration of our map can be written as: 

On = 6o + nro + 0{ne), 

n—1 n—1 

Tn = ro + e'^uJk[v{9k,rk)+ev2{Ok,rk)]+e‘^'^E2{Ok,rk) + 0 {ne^^'"), 
k=0 k=0 

(38) 

where V 2 { 0 ,r) is a given function which can be written explicitly in terms of 
v{9,r) and Si{9,r). 

Recall that IjS is a totally irrational segment ii p/q G I 13 , then |g| > 
where 0 < 2b < (3. 

We recall that we define b = {13 — p) 12 for a certain Q < p < (3. In 
the following we shall assume that p satisfies an extra condition, which will 
ensure that certain inequalities are satisfied. This inequalities involve the 
degree of differentiability of certain functions. We assume that I > 12. 
Then we have that: > 0, lim/_,.oo = 1. Thus, there exists a constant 

R > 0 such that: 

> R > 0, for all I > 12. (39) 

Given (3, satisfying: 

0 < /d < 1/5, (40) 

then we will take p satisfying: 

0 < p < R13. (41) 
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Lemma 5.1. Let g be aC^ function, I > 12. Suppose r* satisfies the following 
condition if for some rational p/q we have \r* —p/q\ < efi then 1^1 > . 

Then for any A such that 


2(3 < A < {I — 1)6 — /3, Q < T = 1 — 2 A < minIA — 2/3, (/ — 1)6 — A — (3, (3} 

and e small enough there is N < such that for some K independent of 
e and any 9* we have: 

„1 N-l 


N / g{e,r*)de+ 

k=0 

In particular, one can choose any 0 < /3 < 1/5, A 
(3/3, b = fi/3. 




= 7/3/3, r 


A-2(3 = 


Proof. Denote go{r) = g(9,r)d9. Expand g{9,r) in its Fonrier series, i.e.: 


9(«.r)=9o(r)+ 

m£Z\{0} 

for some gmif) : M —)■ C. Then we have: 


N-l 


N-l 


J2W + kr’,T‘)-g„(r-))^Y. E 9™(r*)e“”<"‘ 

k=0 m£Z\{0} 


k=0 

N-l 


N 


k=0 l<|m|<[£-'>] k=0 |m|>[£-*'] 

N-l N-l 

9 m{r*)e 

l<|m|<[£-'>] fc=0 


\ „2'Kim6* \ ^ 2nimkr* , \ ^ \ ^ f *\ 2nim(d* +kr*) 

/ . ^ . 9 m\> je 

k=0 


9m{r*)( 


flTTimO* 


- 1 


o27Timr* 


rr + SZ 9m{ry 


l<|m|<[£“*’] 


k=0 |m|>[£-'>] 


PTrim(8*+kr*) 


(42) 


To bonnd the hrst snm in (42) we distingnish into the following cases: 


• If r* is rational p/q, we know that |g| > e 

— If |g| < e~^, then pick N = |g| and the hrst snm vanishes. 

- If \q\ > e then by dehnition of r* for any s/m with \m\ < e ^ 

we have or |mr* — s| > By the pigeon hole principle there exist 

integers 0 < iV = g < and p snch that |gr* — p\ < 2e^. 
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• If r* is irrational, consider a continuons fraction expansion Pn/<in —t r* 
as n —)■ cx). Choose p'/q' = Pn/qn with n snch that g„+i > This 
implies that \q'r* — p'\ < l/\qn+i\ < The same argnment as above 
shows that for any \m\ < we have |mr* — s| > . 

Let N be as above. Then 


9m{r*) 

l<|m|<[£-^] 


^2'KiNmr* 


^27Tim6* 


o2'Kimr* 


- 1 


< 2e^-^ Y. 


We point ont that since g{0,r) is C*, then its Fonrier coefficients satisfy 
\gm{f'*)\ < m 7 ^ 0. Thns we can bonnd the first snm in (42) by: 


Y 3m{r*) e 


^2'KimQ 


* 6 


2'KiNmr* 


- 1 


o2'Kimr* 


(43) 




A-f 


l<|m|<[e“'’] 




(44) 


To bonnd the second snm we nse again the bonnd for the Fonrier coeffi¬ 
cients gm{r*): 


N 


Y 

k=0 |m|>[£-'>] 


<N Y 

|m|>[e“*'] 


(45) 




Clearly, taking r = 1 — 2^4 < minjA — 2f3, {I — 1)6 — A — /3,/3}, and 

□ 


snbstitnting (43) and (45) in (42) we obtain the claim of the lemma. 


Lemma 5.2. Let (3 satisfy (40), and b = {jd — p)/2 with p satisfying (41). Let 


Up be an exit time of the process {0n,rn) defined by (38) from some bounded 
domain ly. Let 6 > 0 be small enough. Suppose that ny > e“2(i-/3)-i-(5^ 
all I > 8 the following holds: 


1. Given two functions h : 
constant d > 0 such that: 


-)■ 


and 5 ^ : T X 


-)■ 


there exists a 


ri/s-l 

Y (^~^"^^h{rk){g{9k,rk) - go{rk)) = 0{e^), 

k=0 

where go{r) = jY(9,r)d9. 
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2. If np < £-2 (i-/3 )+(5^ then given a function : T x M —)■ M and a 
collection of functions : M —)■ M, with ||/?.fc||co < M and ||/ifc+i — 
hk\\c° < for all k, there exists a constant d > 0 such that 


np-l 


hk{rk){g{9k,rk) - go{rk)) = 


k=0 


where go{r) = g{e,r)de. 


Proof. We shall prove both claims using Lemma 5.1 To that aim hx 2/3 < 


A < min{(Z — 1)6 — ,3, (1 — /5)/2}. We note that ( |40| ) and ( |4l| ) ensure that 
2/3 < min{(Z — 1)6 — /3, (1 — /3)/2}, so there always exists such A. 

Since we have ri/j > and we have taken A < (1 —/3)/2 < l — (3 — 5 


(the second inequality being sati shed because 6 is small), we have e ^ < ng. 


Choose N < e ^ from Lemma 5.1 and write np = PgN + 
integers Pg and 0 < Qg < N. Then: 

np-l 

e~^‘'^^h{rk){g{9k,rk) - go{rk)) 

k=0 

Pp-^N-l 

h{rkN+j)ig{9kN+j , VkN+j) - go{rkN+j)) 

k=0 j=0 
3b-1 


for some 


< e" 


(46) 


+ 


j=0 


—Ae 


'(^^-^+j)h(rp^jY+j}(g(9p^jy+j,rp^jv+j) - go(rp^N+j)) 


Let us prove item 1. We shall bound the two terms in the right hand side 


of (46) in a different way. Recall that in the normal form ([1^ we have 

L 


I-)- 


9 + r + eEu{9, r) + euu{9, r) + C>/(£^+“) 
r + e^E2{9,r) + Oiie^+^). 


On the one hand we have that for all k < Pjp, and all j < N: 

r kN+j = rkN + 0{Ne‘^), 

and: 


9kN+j — 9kN + jr kN + 0{N^e). 


Hence: 


(47) 

(48) 

(49) 


e ^’‘^~^^'^h{rkN+j)igi9kN+jAkN+j) ~ go{rkN+j)) 

_ g '^^h{rkN){g{9kN + j^kNAkN) — go{rkN)) + 0{e 
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Then; 


f/3-l N-1 

^^^+^'>h{rkN+j){g{OkN+j,TkN+j) — goi'^kN+j)) 

k=0 j=0 


< e" 


Pfl-i 


N-l 


y~] e ^^hjrkN) 'y~^ (gjOkN + jrkN, tun) — go{rkN)) 


k=0 


j=0 


Pb-1 


+KN^e^ 


-Xe^kN 


k=0 


Thus, using Lemma 5.1 we obtain: 

P/3-1 


N-l 


"Y. ® ^^h{rkN) 'Y, jgj^kN + jPkN, TkN) — goi^kN)) 

k=0 j=0 

Pp-l 

< Y e~^"^^^\h{rkN)\ < 


k=0 


for some constants K, K > 0. Moreover, we have: 


Pfl-i 


KN^e^ Y 


k=0 


Thus: 


P/3-1 N-l 

EE^-"‘ i^^+i)}l{rkN+j){g{dkN+ji PkN+j) — go{PkN+j)) 
k=0 j=0 


(50) 


On the other hand we have: 

3-1 


^ g Ae (^/i^+J)h(^rp^jY+j)(g(dp^N+j,rpf}N+j) - goirPfjN+j)) 
j=0 


<e^K sup \h{r){g{e,r) - gQ{r))\Qp < Ke^ 


{e,r)&Ip 


(51) 


In conclusion, using (|50|) and (|51|) in equation (46) we obtain: 


np-l 


^''^’'h{rk){g{9k,rk) - go{rk)) 


k=0 
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Denoting d = niin{2 — A, 1+/3 — 2A, t + (3}, and letting 1 — 2A = t > (3/3. 
the proof of item 1 is finished. 

Now let us prove item 2. The proof is very similar to item 1. We can use 
the same formula (46) substituting e~^^^^h{rk) by hk{rk). Again, using (48) 


and (49) we can write: 


hkN+j{jkN+j){g{dkN+jAkN+j) — goiAkN+j)) 

= hkN(j’kN){g{dkN + jrkN, ^kN) — goi^kN)) + 0{N'^e). 


Then; 


< 


f/3-1 N-1 

EE hkN+j{rkN+j){g{dkN+jAkN+j) — go{f'kN+j)) 

k=0 j=0 

P/3-1 N-1 

y~] hkNjPkN) + jTfcjV, PkN) “ 5'o(’"fcAr)) 


fc =0 


i=o 


+ KN^e^P^. 


Thus, using Lemma 5.1 we obtain: 


Pfl-i 


N-1 


hkN^PkN) '^^{g{0kN + jrkNyPkN) “ 5'o(?"A:v)) 


fc =0 


j=0 


< 


(52) 


where we have used that P/s < njs < £-2+2/3-5 < £“ 2 +^_ same reason 

we have .Thus: 

Pp-l N-1 

EE hkN+j{PkN+j){g{0kN+j,PkN+j) ~ go{PkN+j))\ E 

A:=0 j=0 

On the other hand we have: 

Q/3-1 

hp,3N+j{rppN+j){g{0ppN+j,rPfiN+j) - go{rp^N+j)) 

j=0 


(53) 


In conclusion, using (52) and (53) in equation (46) we obtain: 

np-l 


e 


e ^''^^hk{rk){g{6k,rk) - giiirk)) 


k=Q 




Choosing r = 1 — 2A > 0 the last two terms are the same. In particular, 
A < 1/2 and 2 — A > 3/2. Therefore, the first term is negligible. □ 
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Let ro belong to the TI case. Consider an interval Ijs = {{6, r) E T x M. : 
k ~ ’"ol ^ for some 0 < /3 < 1/5. Denote G N the exit time from Ip, 
that is the hrst number such that {6np,rnp) ^ Ip- 


Lemma 5.3. Let /? satisfy (40), and b = {ft — p)/2 with p satisfying (41). 
Take / : M —)■ M &e any function with I > 12. Then there exists d > 0 such 
that for all X > 0 one has: 


E 


e" Y, 

k=0 


Kfirk) - ( l>{rk)f'{rk) + f'‘f ''{rk) 


-f(r„) = 0(^-1’*“), 


where for E 2 { 0 ,r), defined in (flip, we have 


b{r) = / E2{0,r)d9, cr‘^{r) = / v‘^{9,r)d9. 


Proof. Let us denote: 


np-l 

k=0 

(54) 

First of all we shall use the law of total expectation. Fix a small enough 
5 > 0. Then we have: 



+ E (r; I < e-2(l-/3)+5) < ^-2(l-/3)+5| 

+ E (r; I > £-2(1-/3)-5) > ^-2(l-/3)-5|_ 

By Lemma |B.2| for e sufficiently small and c > 0 independent of £ we have 

P{n;3 < £-2(1-^)+^} < exp . (55) 

Now we write: 


np-l 

e-^^"^^f{rnp) = /(ro) + - e-"^^'=/(r,)) . 

k=0 


Doing the Taylor expansion in each term inside the sum we get: 


np-l 


= /(i-o)+ [-AA-"'’V(rO + e-''“Yh)(ri+,-rO 


k=0 


+ 0{e-^^ 
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Substituting this in (54) we get: 


np-l 


k=0 


V = f{ro) + 

np-l 

E 




/'(rfc)(rfc+i - Tk) + ^^'^f"{rk){rk+i - r^f 




fc =0 


h{rk)nrk) + (r,) 


-1 


+ C>(e-^^V).(56) 


k=0 


We note that using (38) we can write: 


rfc+i - rfc = £i:nfc[n(6'fc,rfc) + £^2(6'^, r^)] + £^^2(6'^, r^) + C>(£^+“), 


and also: 


(rfc+i - Tk) = £ n (4, Tk) + 0(e:'^). 


Thus we can rewrite (56) as: 

np-l 

V = fi^o) + ^ X] e~^^''^f{rk)uk [viOk, r^) + ev2{0k, r^)] 

fe =0 

71,3-1 

[^ 2 ( 6 'fc,rfc) - 6 (rfc)] 

fc =0 
2 ’T’/S-l 

+ |- Y1 [v^(0fc,rfc) -a^(rfc)] 

fc =0 
np-l 

k=0 


(57) 


Now we distinguish between the case 5 “ 2 (i-/ 3)+(5 < 77 ,^ < £;-2(i-/3)-5 

Consider the former case. First, we show that the last term 
in (57) is for some d > 0. Indeed, 


71,3-1 

C>(e-A£2fc^2+a^ 

k=0 


< Ke^+^n^ < Ke‘^^+^, (58) 

where d = a — 6 > 0 due to smallness of 6 , and K is some positive constant. 
Now we use item 2 of Lemma 5.2 in (57) twice. First we take hk{r) = 
g-Ae and g{ 6 ,r) = E 2 { 6 ,r), and after we take hk{r) = e~^‘^ ^f"{T) and 

g{0,r) = v‘^{9,r). Then, recalling also (58), equation (57) for £- 2 (i-/ 3)+<5 ^ 

< e:“2(i-/3)-<5 yigig^g- 
np-l 

V = f{ro) [v{0k, Tk) + ev2{9k, r^)] + C>(£^^+^). (59) 

k=0 
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Now we focus on the case rip > e The last term in (57) can be 

bounded by: 


np 


-1 


Eo( 


g— Ae^fc^2+a^ 


k=0 


np 


-1 


< Ke 


2+a 


E 

k=0 


—\e^k 


= 


X _ g-Ae^n^ 

1 — 


< Kxe% 


(60) 

for some positive constants K and Kx. Similarly to the previous case, we 
use item 1 of Lemma 5.2 in (57) twice. First we take h{r) = f'{r) and 
g{0,r) = E 2 { 0 ,r), and after we take h{r) = f”{r) and g{0,r) = v‘^{6,r). 
Using this and bound (60) in equation ( [5^ , we obtain the following bound 
for np > e:“ 2 (i-/ 3 )-< 5 . 


np-l 

V = /(^o) + ^ X] e~^^^’'f{rk)ujk [v{9k, Vk) + ev2{9k, r^)] + (61) 

k=0 

Now we just need to note that since cuk is independent of and 9k, we 
have for all A; G N: 


E{ujkf{rk)[v{9k,rk) + ev2{9k,rk)]) = 

E{uk)E{f{rk)[v{9k,rk)+ev2i9k,rk)]) = 0 , 

because ]E(a;fc) = 0. Thus, denoting if we take expectations 

in (j5^ and (61) and use the total expectation formula, it is clear that: 

Hv) - f{ro) 

^ E I ^ e~^^^^f{rk)ujk [v{9k, Vk) + ev2{9k, rk]) ) E{np = n} 


= £ 


nGN 

n>n£ 


,k=0 


+ O (g2/3+rf^p|g-2(i-/3)+5 ^np< e-sii-/?)--?} p 0{e‘^)F{np > £-2(i-/3)-'5} 


'n-l 


= e 


^ ^ e ^^^''E(/'(rfc)a;fe [v{9k, rk) + ev2{9k, Vk)]) F{np = n} (62) 


nGN \ k=0 
n>ne 


+ O (g2/3+rf)p|g-2(l-/3)+5 < < ^-2(l-/3)-5| ^ 0{e'^)F{np > £-2(1-/3)-5| 

= O (g2/3+d)p|g-2(l-/3)+5 ^np< £-2(1-/3)-'5} p 0{e^)F{np > e-2(l-/3)-'5} 


By Lemma B.2 there exists a constant a > 0 such that: 

r{n0 > = o (e^A j , 


(63) 


Clearly, if (63) is true then P{n^ > e: '^} is smaller than any order of 

£ and then (62) finishes the proof of the lemma. 


36 















To prove (63), let us define ns := Define also n, := 

z = 0,..., Clearly, if > £-2{i-/3)-<5^ kn^+i ~ ’"nj < 2^^ for all i. 

In other words, we have that: 


P{n^ > £ < P{|r„.^i - r„J < 2e^ for alH = 0 ,..., [e '^]} 


[£-«] 

JJ - r,, 


m+1 ' n^ 


< 2 eP. 


(64) 


i=0 


where in the last equality we have used that — r„. and are 

independent if i 7 ^ j. 

Now, take any i. Then: 


n-i+i-l 


Tui+i -rni= e ^ ujkv{9k, Vk) + 0{e^{ni - ni+i)). 


Note that e^{ni — nj+i) = e^[e . Thus: 

»4i + l-l 


uJkv{dk,rk) 


- 0{e^^) < -r„,|. 


(65) 


As a consequence, if - rn,\ < 2e^ then £ |X)r=ti) ^^kv{6k,rk)\ < 3e^. 
Indeed, if this latter inequality does not hold, then: 


V ^kv{0k,rk) 


- 0(e^l>) > Ss'’)! - 0(sl>)) > 2s'’ > |r„.^. -r„,|, 


which is a contraditciton with (65). In other words: 


r-n,+i -rn^< 2e^} < P <j £ 


rii+i-l 

ujkv{9k,rk) 

i=ni 


< 3e^ 


Now by Lemma B.2 that (uj+i — Ui) X)fc=n ^kv{9k, Vk) converges in dis¬ 


tribution to ^ ~ A/ ( 0 , c^) for some c > 0 . 

Thus, using that rij+i — Ui = as £ —)■ 0 we obtain: 


P < e 


fii+i-i 

y^ ujkv{9k,rk) 

i=ni 


<3£n=P{|e|<3} + o(l)<p<l, 


for some constant p > 0. Then, using this in (64) we get: 

-5 


> £-2(1-/5)-5| < pl/e^ 


Dehning a = — logp > 0 (because p < 1) we obtain claim (63). 


□ 
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5.2 The Imaginary Rational case 

In this section we deal with the imaginary rational case. The ideas are 
basically the same as in the TI case. Recall that after performing the change 
to normal form the n-th iteration of onr map can be written as: 


On = 6Q + nrQ + 0{ne), 

n— 1 

ro + e'^ujk[v{9k,rk) +ev2{0k,rk)] 
Efc=(! E2i9k, Tk) + C>(ne2+“), 


r„. = 


( 66 ) 


where V2{9,r) is a given fnnction which can be written explicitly in terms of 
v{9,r) and Si{9,r). 

We also recall that given an imaginary rational strip Ip there exists a 
nniqne r* G Ip, with r* = p/q and 1^1 < Moreover, for all tq G Ip we 
have |ro — r*| < e^. Then by ( 66 ) for any n < np we have: 


r„ = 


Dehne 


6*0 + nr* + C>(ne^), 

(67) 

r* + 0{e^). 

1 

r) = - ''^^g{9 + ir, r). 

^ i=0 

( 68 ) 


Lemma 5.4. Let jl satisfy (40), and b = {/3 — p)/2 with p satisfying (41). Let 
np be an exit time of the process {9n,rn) defined by (38) from some hounded 
domain Ip. For all I > 1 the following holds: 


1. Given two G functions h : 
constant d > 0 such that: 


-)■ 


and : T X M —)■ M, there exists a 


np-l 

^ e~^''^^h{rk){g{9k,rk) - go{9k,rk)) = 

k=0 


2. If np < then given a G function : T x M —)■ M and a collection 

of functions : M —)■ M, with ||hfc||co < M and \\hk+i — hk\\co < 

for all k, there exists a constant d > 0 such that: 

np-l 

^ hk{rk){g{9k,rk) - go{9k,rk)) = 

k=0 
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Proof. Let us start with item 1. Write = Ppq + for some integers Pp 
and 0 < < g. Then: 




71/3-1 

Y e~^‘''^k(rk)(g(0k,Tk) - go{0k,Tk)) 

k=0 


< 

Pp-1 q-l 

EE 

k=0 j=0 

e ^ h(^rkq+j){g {Okq+j, Tkq+j) g^iOkq+j iTkq+j)) 


(69) 


+ 


Q ^-1 

j=0 


g ( P'i^Pli{rPi^qj^j)(^g{6ppqj^j.irppq-\.j) go(^P;3q+j5 ^P/3(}+j)) 


On the one hand, we note that by (|67|) and j < q < e ^ we have 


dkq+j — dkq + + 0{£^ ^), 


Tkq+j =r* + 0{£^). 

Then for all k < Pp\ 

e ^ h{rkq+j)(^g{Okq+jiT'kq+j) 9oif^kq+j I'^kq+j)') 
= e-^^"^^h{rkq){g{ekq+3r\r*) - goiOkqP jr^r*)) + 


Then: 


Pp-^ q-l 

EE'- ^ ^^h(rkqpj){g{Okq+j, Tkq+j) gQ^Okq+j, Tkq+j)) 
k=0 j=0 


< e" 


Pr-I 


q-l 


e ^""'“^hirkN) Y^siOkq + jT*,T*) - go{0kq + jT*,r*)) 


k=0 


j=0 


Pr-1 


+Ke^-^^ Y 


^—Xe^kq 


(70) 


k=0 


Now, recalling that r* = p/q, by the dehnition (68) of go{0,r) for all k < Pp 
we have: 


q-l 

Y^3{0kq+ jT*,r*) 


j=0 


q-l 

go{0kq+ jT*,r*)) = Ygi^kq,T*) - qgo{0kq,T*) 


j=0 


0. 


Moreover: 


Pp-i 

^^3-26 ^-Xe^kq ^ ^1-26 

k=0 
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Using these estimates (70) yields: 


Pfi-l q-l 

^ ^ h{T}^q^j^{^g{6i^q^j ^T}^qj^j') Qoifikq+jiPkq+j)')^ ^ 

k=0 j=0 

Ke^-'^K 

Note that 1 — 26 = 1 — /3 + p>0, since (3 <1 and /3 — 26 = p > 0. 

On the other hand, we have: 

Q/3-1 

^^P^^^^h{rppN+j){,g{,PpN+j,rppN+j) - go{rppN+j)) 

j=0 

<e‘^K sup \h{r){g{9,r) - go{9,r))\Q^ < Ke^-\ 

<,e,r)Elp 


(71) 


(72) 


Clearly, 2 — 6 > 0. Substituting ( |71[ ) and ( |72[ ) in ( |69[ ) yields item 1 of the 
Lemma. 

Now let us consider item 2. If we take equation (70) and substitute 




< 


h{rk) by hfc(rfc), we can write for all k < Pf 

f/3-l q-l 

EE hkq+j{Pkq+j){g{^kq+j-i ^kq+j) goi^^kq+ji ^kq+j)) 

k=0 j=0 

P/3-l q-l 

^ hkq(rkq) ^(g(% + jr*,r*) - go{9kq + jr*,r*)) 

j=0 


(73) 


k=0 


+ Ke^-^’^Pg. 


q-l 


Again: 

^{g{0kq+ jr*,r*) - go{9kq + jr*,r*)) = 0, 

j=0 

and since Pg < Uf^ < £-2{i-/3)-<5 ^nd 6 = {(3 — p)/2 

< A:£^+2/3-2b < 

Then we have: 

Pfl-l q-l 

EE^- ^ '^^h(r/jg_|_j)(p( 0 fcij_i_j,rfcg_i_j) po(^fcij+j)’^fcg+j)) 

k=0 j=0 

On the other hand we have: 

Q/9-1 

hppN+j{r p^Ar+j)(p(p3Ar+j, Tp^N+j) - go{rPpN+j)) 


< 


(74) 


j=0 
,2 


< e'KQf^ < Ke^-\ 


(75) 


40 



















We note that 2 — b — 2f3 > 0 for 6 = (/3 —p)/2 and /3 < 4/5 and that 1 —/3 > 0 
if < 1. Thus, taking /? < 4/5 bounds (74) and ( [75| ) prove item 2 of the 
Lemma with d = min{2 — b — 2/?, 1 — /3, p} > 0. □ 

Let To belong to the IR case. Consider an interval J /3 = {( 6 ', r) G T x M : 
k ~ ’"ol < for some 0 < /? < 4/5. Denote G N the exit time from Ip, 
that is the hrst number such that {6np,rnp) ^ Ip- 

Lemma 5.5. Let / : M —>■ M &e any function with I > 3. Then, there 
exists d > 0 such that for all X > 0 one has: 




71,3-1 

E 




e 2 _^e 

fc =0 
-/(ro) = 


- ( &( 6 'fc,rfc)/'(rfc) + 


(r^{0k,rk) 


firk) 


where: 


*7 1 1 ^ ^ 

b{9, r) = - ^ 7^2 (d + ir, r), 0-^(9, r) = - ^ v^(9 + ir, r). 


q 


q 


2=0 ^ 2=0 
Proof. Let us £x any 0 < d < 1/6. Again, denoting: 




71,3-1 

^ e 

k=0 


— Xe^k 


><f{rk) - { b{9k,rk)f'{rk) + 


(J^{9k,rk) 


/"(rk) 


(76) 


and using the law of total expectation, we have: 

E (p) = E (^r? I £-(1-/5)+-^ <np< e-2(l-4)-<5^ p|^-2(l-/3)+5 < < £-2(l-/3)-5| 


+ E (^ri\np < e 
+ E(^y\np> e-2(i-/3)-5^ 


< e 


- 2 ( 1 - 


/3)+<5j 


(77) 


> e 


-2(l-/3)-5 


}• 


As in the proof of Lemma |5.3[ for e: sufficiently small by Lemma B.2 and 
some C > 0 independent of £ we have 


< < exp (-T) , 


and, thus, (77) yields: 


E (p) = E(p I <np< e-2(l-/3)-5)p|^-(l-/3)+<S <np< £-2(l-/3)-<5| 

+E (p I > £-2(1-/3)-5) > ^-2(l-/3)-5|_ 
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Now, we can write: 


np-l 

= /(ro) + f{r,) + f {r,){n+, - n) 


k =0 


+ ^e-"^^V"(rfc)(r,+i - r,)^ + 


and then ( 17 ^ can be rewritten as: 


V = firo) + 

71/3-1 

E 

k =0 


e ^'’V'(ft)(’')d+i - ri) + h *'’‘/"(ri)(r/i+i - rt)‘ 


J( 79 ) 


np 


-1 




Ae^fc 


k =0 


bi^k,rk)f{rk) + 


o-‘^{6k,rk) 


fin) 




-1 




fe =0 


Now, using (66) we have: 


rfc+i - rfc = + £^2(6'*,, r^)] + £^^2(6'^, r^) + C>(£^+“), 


and: 


(rfc+i - rfc)^ = e^v^Ok, Vk) + 0(e:'^). 


Thus, ( 79 ) writes out as: 


n/s-l 

V = f{ro) + ^ e~^^^’'f{rk)eu}k [v{9k, Vk) + ev2{0k, r^)] 

k=0 

77/3-1 

[^ 2 ( 6 'fc,rfc) - 

fc =0 
2 77/3-1 

k=0 

77/3-1 

+ ^ C>(e-^"'V+“). 

k=0 


( 80 ) 


Now we distinguish between the cases e“2(i-/3)+5 < < £“2(i-/3)-<5 

> £-2(i-/ 3)-<5_ shall start assuming that £“2(i-/3)+<5 < 77,^ < £;- 2 (i-/ 3 )- 5 _ 
As in the proof of Lemma 5.3 we have: 


71/3-1 

^ C>(e-^"'V+“) 

k =0 


< < A:£2/3+a-5 


( 81 ) 
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We note that 1/6 — 6 > 0 since we have taken 6 < 1/6. Now we nse 
item 2 of Lemma 5.4 in (80) twice, taking hrst hfc(r) = and 

g{ 6 , r) = E 2 { 6 , r), and after hk{r) = and g{6, r) = v‘^{6, r). Then, 

nsing also (81), eqnation (80) yields: 

np-l 

h = /(^o) + 'Yh e~^^^’'firk)eujk [v{9k, r^) + ev2{9k, r^)] + (82) 

fc =0 


for some snitable d > 0. 

Now tnrn to the case rijs > £“2(i-/3)-5^ 
bounded by: 


The last term in (80) can be 


np 


-1 


Eo(< 


—X€^k2-\-a 


^i+a) 


k=0 


np 


-1 


< ^ = Ke^ 


k=0 


X _ g-Ae^n^ 

1 — 




(83) 

for some positive constants K and K\. Then, using item 1 of LemmajT^twice 
(hrst with h{r) = f'{r) and g{9,r) = E2{9,r), and later with h{r) = f"{r) 
and g{9,r) = v'^{9,r)), we obtain the following bound for np > £“2(i-/3)-(S. 

np-l 

V = fi^o) [v{9k, Vk) + ev2{9k, r^)] + C>(e:'^). (84) 

k=0 

To hnish the proof, we follow the same steps as in the proof of Lemma 
TSjand ([78| yields: 


E(??) - /(ro) 


'n —1 


= f\rk)eu;k [v{9k,rk) +ev2{9k,rk]) | ¥{np = n} 

ngN \k=0 J 

+ C)(£2/3+'i)p{g-(i-/3)+'5 <np< g-sii-/?)--?} + 0{e'^)'¥{np < 

= C)(£24+'i)^ 

where by Lemma |B.2 for some C* > 0 independent of e 


> e 


-2(l-/3)-5x _ 


c 


} = exp ( ) , 


so that it is smaller than any power of e. 


□ 


5.3 The Real Rational case 

Here we study the system in the RR case, in the subdomain: 

\r — p/q\ < . 
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Figure 3: Level sets of the pendulum pjq = 1/3. 


From the Normal Form Theorem, in the Real Rational strips the system 
takes the following form: 

0i= Oo + ro + e [EM(6'o,p/g) - Eu(6'o,p/g) + Eup,q(6'o,p/g) + EsiOo)] 

+ eujou{eo,p/q) + 0{e^^^), 
ri = ro + eEvp^g(0o, ro) + £CJov(0o, ro) 

+ £^'^^a;oU2(6'o, ro) + e^^^B4(0Q, Tq) + 0(£^~^'^), 

where V2{0,r) can be written explicitly in terms of v{d,r) and Si{6,r). The 
function E 4 is such that ||F^4||co < K. We point out that it is a rescaled 
version of the function F^4 appearing in the Normal Form Theorem. 

Recall also that: 


max{||Eu||co, ||Eu||co, ||Eup_g||co, ||.F3||co, \\u\\co, ||u||co} < K. 


Moreover, we have dehned V 2 also in such a way that ||u2||co < K. 
First, we switch to the resonant variable: 


f = r — p/q, 0 < |f| < 


With this new variable, system (85) writes out as: 


01 = 00+Plq + e [Eu{0o,p/q) - Ev{9o,p/q) +Evp^q{0o,p/q) + Es{6o)] 

+ fo + ea;ou(6»o,p/g)+ 0(e^+“), 
fi = fo + eEvp^q{0o, fo) + e / £^4(6*0, fo) + ea;or(6»o, fo) 

+ £^/'^uoV2i9o, fo) + C>(e^+“), 

where: 

'D(6'o,fo) = v{9o,fo+p/q), V 2 { 0 o,ro) = V2{9o,fo+p/q), 

Evp^q{9o,fo) = Evp^q{9o,fo + p/q), £4(6*0 ,fo) = £4(6*0, fo + p/g), 


44 



From now on, we will abuse notation and drop all hats. We are interested in 
the g—th iteration of map (86), which is given by: 


Oq = 0o + qro+ 

g-l 

+ e'^[Eu{ek,p/q) - Ev{ek,p/q) + {q - k)Evp^q{ek,p/q) + ^3(^fc)] 

k=0 

+ [u{0k,p/q) +v{ek,p/q)] + C>(£^+“), 

q-l 

Tq = ro + e'^ [Evp^qiOk, Vk) + Vk)] 

k=0 

+ ^ ELo Tk) + e^/‘^V2{0k, Tk)] + C>(£^+“). 

Taking into account that q is bounded for 0 > i < g we have: 

9i = So + i(plq + r„) + 0(e), r, = rj + 0(e). 


(87) 


we can rewrite the last system as: 

Oq= 00 + gro + eEu^'i\eo) + eu^‘^\eo,uj^o) + 

Tq = ro + £En('?)(6»o,ro,£) + ro, eng) + i^o) + 0{e‘^). 

where u^l = {uqk, ■ ■. ,u}qk+q-i) and: 

q-l 

Eu^’i\d) = E [E-u(6' + ip/q,p/q) — Ev{9 + ip/q,p/q) + 

i=0 

+ (g - i)Evp^q{9 + ip/q, p/q) + E^{9 + ip/q)] , 

q-l 

u^^\9,ujI) = ^{q-i)ujqk+i[u{9+ ip/q,p/q)+v{9+ ip/q,p/q)\, 

i=0 

q-l 

Ev^'^\9,r,e) = '^\Evp^q{9 + i{p/q + r),r) + + i{p/q + r),r)\ , 

i=0 

q-l 


M) 


,(<?) 


{9,r,ujl) = ^ujqk+iv{9+ i{p/q + r),r), 


i=0 

q-l 


{9,r,ujl) = ^ujqk+iV2{9+ i{p/q + r),r). 


i=0 


Introduce a rescaled variable r = R^/e. Then {9i,Ri) are dehned using 
system (87) with the corresponding rescaling. This can be rewritten in the 
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following way, where we just keep the necessary ^-dependent terms: 


Oq = 0^ + qRoV^ + eEu^^\do)+eu^’^\9o,u:l) + 0{e^+^), 

Rq = Ro + e^/‘^Ev^'^\9o,RoV^,0)+e^/h^'^\eo,RoV^,ul) 

+£n2'^^(6'o, 0, cjq) + 0{e^^'^). 

Consider the Hamiltonian: 

r2 1 rd 

H{9,R) = - -/ En('')(s,i?V£,0)ds. (89) 

2 Q Jo 

Let Hfi • ^(9qn:Rqn)' Study the process. (9q<fi^Rn) • Lr(^g^, 

where {9qn, Rqn) is the process obtained iterating ( [^ n times. One can see 
that: 


H, = Ho+ ^eRov^\9^,R,^e,ujl) 

+ ^F{ 9 o, Ro, £) + ^ 0 ( 00 , Rq, cjq , e) + 0 {e^^°‘), 


(90) 


where F and G are: 

F{9,R) 


(0,0,0) Eu^‘i\9) 

Q 

1 

-En('')(0,0,0) / Ev^'^\s,0,0)ds 
Q Jo 

^R^deEv^'^\9,0,0) + ^ (En('')(0, 0, 0))' 


1 


q-l 




G{9,R,ul) = 


i=0 

-En('')(0,0,0) u^'^\9,ul) 

-n('?)(0,O,a;^,O) / a,En('^)(s, 0, 0)ds 

^ Jo 


q-l 


^qk+j v{9 + ip/q,0)v{9 + jp/q,0) 


i,j=0 

i¥=j 


Ev^'^\9,0,0)v^^\9,0,ul) 


q-l 


+ R UJqk+iV2 (0 + ip/ g), 0) • 


i=0 


We note that since \R\ < Ki we have: 

||F||co<iL, ||G||co<iL. 
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Moreover, one has that for all /c > 0: 


^{G{6qki Rqk,^^) — 0 . 


In terms of the variable if, the {|i?| < Ki\ region can be written as: 


Irr := {if e M : |if | < ffg}, 


(91) 


for some constant K^. Denote by n* the stopping time of leaving Irr for the 
first time. 


Lemma 5.6. Let / : M —)■ M &e any G function with I > 3. Then there exists 
d > 2/5 such that for all X > 0 one has: 


n*—l 

k=0 


XfiHk) - biOqk, R,k)fiHk) - 


-/(iio) = 0(7), 


where: 


q-l 

b{e,R)=F{e,R), a\e,R) = Rl,J2v\e + tp/q,0). (92) 

i=0 

Proof. Depending on properties of the averaged potential , Ry/e, 0) 

we decompose the region {|if| < ffi} into several domains as follows. 

In the interval [—ffi,ffi] there are finitely many critical valne of if dne 
to the fact that is a trigonometric polynomial in 9. By assnmption 

[H5] all critical points are nondegenerate. Therefore, there are hnitely many 
of them. Consider a y^-neighborhood of each one of them. 

There are level sets that are not intersecting these neighbonrhoods of the 
critical points. Call these level sets regular. The others levels intersect these 
neighbonrhoods. Call these level sets nearly critical. Notice that the regnlar 
level sets are no necessarily connected (see Fig. |^. 

Each connected component of a family of ovals gives rise to a segment 
and the nnion of snch segments gives rise to a diffnsion on a graph as in 
Freidlin-Wentzell [20], Sect. 8. 

Now we stndy separately three regimes: 

• a y^-near a critical point; 

• iterates with an initial condition on a near critical level set; 

• iterates with an initial condition on a regnlar level set; 
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In the first case note that n* < oo almost surely. Indeed, on the one hand 
if {6*, 0) is a critical point of the dominant part of the deterministic system, 
i.e. 0, 0) = 0,0) = 0, then the drift of the process is given 

by: 

0 ) = 2 

i=0 


On the other hand, if {0,0) is not a critical point, then the drift of the 
/^-component of the process {6qn,Rqn) defined in (88) is precisely, which is 
0,0) 7^ 0. In conclusion, the process does not get stuck at i? = 0. 
Since away from this zone the diffusion coefficient a{9, R) is nonzero, n* must 
be hnite with probability one. 

In the second case denote: 


rj = e-"-7(^^n‘) + 

n*-l 


k=0 

Now we write 


—Xek 


Xf{Hk)-b{9qk,Rqk)f{Hk)- 


^ {9qk, Rqk) 


r{Hk) 


(93) 


n* — l 


= f{Ho) + [-Ae-^^7(^fc) + - Hk) 


k=0 


1 


+-e-^^’^r{Hk){Hk+i - Hkf + 0{e-^^h^R) 


so that (93) writes out as: 
f{Ho)+ 


V = 

n* — l 


E 

k=0 

n*-l 

E 


1 


e-^^’^f{Hk){Hk+i - Hk) + - H^)^ 


(94) 


-£2^ e 
k=0 


— Xek 


b{eqk,Rqk)f'{Hk) + 


^ {9qk, Rqk) 


f"{Hk) 


n* — l 


+ 0(e-^"^e3/2)_ 


A;=0 


Now, using (90) it is clear that: 


Hk+l Hk y/^Rqk^^'^^ {9qk, Rqk's/^, ^k')R 

eF{9qk, Rqk) + £G{9qk, Rqk, ojD + 0 ( 5 ^+“). 


(95) 


Moreover, we have: 

{Hk+i - Hk)^ = ER%{v^^\6qk,Q,u:l)f + 0{e^R) 

q-l 

= £Rlk'^v^{9qk + ip/q.O) + eGo{eqk,Rqk,i^l) + 0{e^R), 
i=0 


(96) 
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where: 


GQiOqki Rqki ^k) 
q-1 q-1 

= E E ujqk+li^qk+jviO + Ip/q, 0)v{e + jp/q, 0). 

1=0 j=l+l 

We note that K{Go(^k, Rqk,^k)) ~ 0- 

Using ( p5| ) and ( [^ , equation (94) writes out as: 

n*-l 


p = f{Ho) + .reY,^~''"''f\Hk)RqkV^^\eqk,V~eRqk,ool) 

k=0 

n*-l 

+e Y, e-^^^f'{Hk) [F{eqk,Rqk) - b{eqk,Rqk)] 


k=0 

n*-l 


+ 5 E 


k=0 

n*-l 


q-1 


R^k Y “ (^‘^{Oqk, Rqk) 


i=0 


«E 

k=0 
n*-l 

+ Eo( 

k=0 


— Xek 


rm, 


f{Hk)G{9qk,Rqk,0jl) + '-Y^Go{9qk,Rqk,Ujl 


^-Xek^l+a^ 


Now, on the one hand by dehnition of b{6, R) and R) it is clear that: 

F{dqk, Rqk) - h{6qk, Rqk) = 0 

(97) 

Rlk + w/R 0 ) - a'^i^dqk, Rqk) = 0 . 

i=0 


On the other hand, the last term in (97) can be bounded by: 

n*-l 


Eo( 


g-Art^l+a) 


k=0 


n*-l 


< Y 


1 — e 


—\sn* 


k=0 


1 — e 


—\e 


< I<xe\ 


(98) 


for some positive constants K and Kx- Using (97), (97) and (98) in equation 
(97) we have: 


n*-l 


p = f{Ho) + i/eY^~^"''f\Hk)RqkV^‘^\eqk,i/^Rqk,uj\ 


(99) 


k=^ 


n*—1 




—\ek 


k=0 


f{Hk)G{eqk,Rqk,ul) + 


nHk 


'Goidqki Rqki ^ 


1) +0{e^). 
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Thus, to finish the proof, we just need to use that: 

E V^Rgk^ooD) =E{Go{e,k,R,k,ujl)) =E{G{e,k,R,k,ujl)) = 0, 

and that Rgk and Hk are independent of ool- Using these facts in ( p^ , one 
obtains straightforwardly: 

E(r/)-/(Tro) = 0(0, 

and the proof is finished. □ 


5.4 Transition Zones 

Here we study the system in the RR case, in the subdomain: 

< |r — p/g| < K 2 e^^^, 

for certain constants Ki and K 2 . From the Normal Form Theorem, in the 
Real Rational strips the system takes the following form: 


di — 9Q + rQ + 0{£), 

ri ^ tq + eEvp^g{0o,ro) + eujov{0o,ro) + 


( 100 ) 


Recall that: 


||Enp,g||co < K, ||n||co < K. 

For our purposes, it will be more convenient to work with the variable: 
f = r — p/q, < |f I < K 2 E^I^. 

With this new variable, system (100[) writes out as: 


Q\ — 9q + p/q + f-Q-\-0{£)^ 

fi ^ fo + eEvp^giOo, h) + eujov{0o, ro) + 


( 101 ) 


where: 

v{9o,fo) = v{9o,ro+p/q), Evp^q{9o,fo) = Evp^q{9o,ro + p/q)- 

From now on, we will abuse notation and drop all hats. We are interested in 
the g—th iteration of map (101), which is given by: 

9q = 9Q + qrQ + 0{£), 


Tq = 


q-1 q-1 

ro + £y^Evp^q{9k,rk) + e'^Ukv{9k,rk) + 

k=0 k=0 
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Note that we have used that |g| is bounded. Moreover, taking into account 
that for i < q, we have: 

= 00 + i{p/q + ro) + 0{e), n = vq + 0(e), 

so that we can rewrite the last system as: 

0g = 00 + 

( 102 ) 

Tq = ro+ £En('^)(0o,ro)+£n^'')(0o,ro,a;g) + 

where we use the notation = {ojkq) ■ ■ ■ )^kq+q-i) and: 

q-l 

En^'^^(0,r) = y~~^Enp^q(0 + i{p/q + r),r), 

i=0 

q-l 

v^'^\9,r,ul) = y^^Ukq+iv{9+ i{p/q + r),r). 

i=0 

We point out that for any k > 0: 

E {v^‘^\eqk,rqk,ujl)) = 0. 

Now let us consider the following function: 

En*^'^^(s, r)ds. 

In this section we will study the process: 




Hn • H(9qYi^Tqn)^ 


where {9qn, ^ qn) 
easily see that: 


is the process obtained iterating (102) n times. 


One can 


Hi = Ho +erov^'^\9Q,ro,u:l) + 0{erl). (103) 

Now, given a constant 0 < p < 1/6, we want to study the process Hn in 
the following if—strips: 

/p(ro) = {if e M : |ii - fiol < Irole'-^- (104) 

We stress out that the width of these strips depends on the initial condition 
To and £. To avoid this, we can dehne the process: 

Hn := |ro|-Vi+"ff„. 
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The process Hn is defined through: 

H, = ^o + ^V‘')(0o,ro,a;o") + O(£lro|). (105) 

Clearly, if we denote by Up the hrst exit time of the process from the strip 
Jp(ro), it is also the hrst exit time of the process Hn of the strip: 

/ = {if eM : |hf-i7o| < !}• (106) 


Lemma 5.7. Let / : M —)■ M &e any function with I > 3. Then there exists 
d > 0 such that for all X > 2/5 one has: 


rip-l 

k=0 

-fiHo) = 0(£"), 


-Ae^Pfc 


XfiHk) - 


(^‘^{OqkHqk) 


nSk 


where: 


q-l 


a 


'{6, r) = ^ v^{d + i{p/q + r), r). 


i=0 


Proof. Let us denote: 


rip-1 

V = e 


^-Ae^Pr 


2p ^ ^-Ae2pfc 

k=0 


\f{Hk 


(J^iOqkHqk) 


ns,) 


■ (107) 


Now we write: 

.^-Xe^Pr 


e-^^^^-pf{Hnp) = 

f{So) + 


rip-l 




k=0 


+-e-^^""’^nH,){H,+i - sn+o{e-^n^n 


so that ( 107| ) writes out as: 
V = fiSo) + 

rip-l 

E 


(108) 


k=0 




j : e 


-Xe^Pk^'^i^qkHqk) 


ns,) + ^ o{e-^n^n- 


k=0 


k=0 


Now, using (105) it is clear that: 


Hk+i - Hk = ef'v^‘^\9gk,rgk,ujl) + 0{e^\rgk\). 


(109) 
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Moreover, we have: 

{Hk+,-Hk)^ = e^P{v^<^\e,k,rgk,u:l)y + 0{e^^\r,k\) 

q-l 

= g 2 p v‘^{eqk + i{p/q + rqk),rqk) 
i=0 

+e‘^^Go{9qk,rqk,(jj'l) + 0{e^^\r qk\) (HO) 

where: 

q-l q-l 

Go{e,r,ujl) = 2 EE 0 Jkq+i 0 Jkq+jv {6 + l{p/q + r),r)v {6 + j{p/q + r),r). 

1=0 j=l+l 

We note that E(Go(6'gfc, cj^)) = 0. 

Using ( 109| ) and ( |110 ) and noting that e^'^\rqk\ < s'^\rqk\., eqnation (108) 
writes ont as: 


TLo-l 


V = 


im+e'’j2 


-Xs'^Pk p// 


f\Hk)v^‘^\0qk,rqk,UJ 


k=0 


^2p ^ 

k=0 
_2p "'P~^ 


■q-l 


'^v'^idqk + i{p/q + rqk),rqk) - cr‘^{9qk,rqk) 


i=0 


E + Y, 0(e-"'’"‘e'’|r,t|). (Ill) 

/c=0 k=0 

Now, on the one hand by dehnition of a^(9,r) it is clear that: 

q-l 

'^v^{9qk + i{p/q + rqk),rqk) - (T‘^{9qk,rqk) = 0 . ( 112 ) 

i=0 


On the other hand, the last term in (111) can be bonnded by: 


Eo(' 


^—Xs^Pkp 


^^Kkl) 


k=0 


1 — e 




-Xe^Pr. 


k=0 


1 - 
=U 

< Kxe^^^-P, (113) 

for some positive constants K and K\, where we have used that \rqk\ < 
if 0 < gfc < Up. Using (112) and (113) in equation (111), and denoting 
(i=l/6 — p>0, we have: 


rin—l 


c-2p '^P~^ 

^ -Xe^Pkrili 


k=0 


WE 


'S"(Ht)G40,k,r,t,ojl) + 0(6“). (114) 


k=0 
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Thus, to finish the proof, we just need to use that: 


E {v^‘^\6qk,rqk,U}l)) = E {Go{6qk,rqk,(^l)) = 0. 


Indeed, using this fact in (114), one obtains straightforwardly: 


E(<7) - f(S„) = 0{e^), 


and the proof is finished. □ 

A Measure of the domain covered by RR and 
IR intervals 

In this section we show that, with the right choice of 6, the measure of the 
the union of all strips of RR and IR type inside any compact set: 

A /3 = Ufc/^ C T X R strips of width 2e^ 

goes to zero as e —?• 0. 

In fact, we will do the proof for A = [0,1]. The general case is completely 
analogous. Let us consider: 

7^ = {p/q eQ : p<q, gcd(p, q) = 1, q < C [0,1], 

where ^max = and: 

Tlq = {p/q eQ : p < g, gcd(p, q) = 1}. 


Finally we denote: 


I^^ = {I|C[ 0 ^] : 3p/qennlp}. 

Lemma A.l. Let p he fixed, < p < (3, and define h = {(3 — p)/2. Then, 
for each Ip such that there is at most one rational p/q satisfying \q\ < 
the union Iqi has the Lebesgue measure p^In) < ond, therefore, as e ^ 0: 

T{In) -e 0 , 

where p denotes the Lebesgue measure. 

Proof. On the one hand, suppose that p/q & ly, q < s~^. Then, for all 
p'/q' G 1 / 3 , with p' and q' relatively prime and p'/q' ^ p/q, we have: 

> \p/q-p/q\ 

qq' q' 
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Therefore: 


> ^-P+b ^ ^-fe-3p/2 ^ 

SO the first part of the claim is proved. 

On the other hand we note that, if qi ^ q^-, then = 0. Moreover, 

it is clear that <7 — 1 (and if q is prime then = 7 — 1, so that the 
bound is optimal). Therefore we have: 

Qinax ^max 2 

#K < < ^9 - 1 = ^ < £-“ 

q=l q=l 

Since /i(// 3 ) = one has: 

0 < /i(/7^) = 

so that the second claim of the lemma is also clear. □ 

B Sufficient condition for weak convergence 
and auxiliary lemmas 

In order to prove that the r-component exhibits a diffusion process we need 
to adapt several lemmas from Ch. 8 sec. 3 [20] • We recall some terminology 
and notations (see Ch. 1 sec. 1 IZDl for more details). 

In the notations of section 11.61 we have 

Lemma B.l. (see Lm. 3.1, Let M be a metric space, Y a continuous 
mapping M 1 —)■ Y{M), Y{M) being a complete separable metric space. Let 
{Xf, P^) be a family of Markov processes in M; suppose that the process 
Y{X() has continuous trajectories. Let {yt,Py) be a Markov process with 
continuous paths in Y{M) whose infinitesimal operator is A with domain of 
definition Da. Suppose that the space C[0,cxd) of continuous functions on 
[ 0 , cx)) with values in T is taken as the sample space, so that the distribution 
of the process in the space of continuous functions is simply Py. Let T be 
a subset of the space C(Y{M)) such that for measures pi, p 2 on Y{M) the 
equality f /dpi = f fdp ,2 for all f E implies pi = /X 2 - Let D be a subset 
of Da such that for every / G T and A > 0 the equation XF — AF = f has 
a solution F E D. 

Suppose that for every x E M the family of distributions Q% ofY{X^) in 
the space C[0, cxd) corresponding to the probabilities Pf for all e is tight; and 
that for every compact K C Y{M), for every f E D and every A > 0, 

poo 

EJ / exp(-A() l\f{Y(Xt)) - Af(Y(Xn)] dt ^ /(K(x)) 

Jo 

as e —)■ 0 uniformly in x E Y~^{K). 

Then Q% converges weakly as e ^ 0 to the probability measure Py{x)- 
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In our case Y (M) is the real line. We use a discrete version of this lemma 
in our proof. 

Similarly, to Lemma 3.2 [20] one can show that the family of distributions 
(those of Y (X^) with respect to the probability measures in the space 
(^[OjOo)) with small nonzero £ is tight. Indeed, in our case speed of change 
of / is bounded. Denote H{X) = H{r,9) = r^/2. Then 

• for every T > 0 and 5 > 0 there exists Hq such that 

PUmax \H{X^)\ > Ho} < 6. 

• for every compact subset K (Z A and for every sufficiently small p > 0 
there exists a constant Ap such that for every a E K there exists a 
function /“(p) on D(A) such that /“(a) = l,/“(p) = 0 for p{y,a) > 
P, 0 < fp{y) < 1 everywhere, and fp(Y {X})) + Apt is a sub martingale 
for all £ (see Stroock and Varadhan 1311). 

In the proof we need an auxiliary lemmas. We study the random sums 

n 

Sn = n > 1, (115) 

k=l 

where {a;fc}fc>i is a sequence of independent random variables with equal ±1 
with equal probability 1/2 each and {vk}k>i is a sequence such that 

lim ^ = a. 

n^oo fl 


Here is a standard 

Lemma B.2. {Sn/rA^'^}n>i converges in distribution to the normal dirtri- 
hution A/'(0, (T^). 

Recall that a characteristic function of a random variable X is a function 
^ C given by (j)x(t) = lEexp(RX). Notice that it satishes the 
following two properties: 

• If X, X are independent random variables, then (px+y = ‘Px ■ Py- 

• Pax{t) = Px{at). 

A sufficient condition to prove convergence in distribution is as follows. 

Theorem B.3 (Continuity theorem [6]). Let {X„}„>i,X be random vari¬ 
ables. If {px„it)}n>i converges to pyit) for every t E M., then {X„}„>i 
converges in distribution to Y. 
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A direct calculation shows that 


lim log (j)s„/^{t) = 


2a2 


for all t G M. 


This way of proof was communicated to the authors by Yuri Lima. 
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